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a 
ABSTRACT 

Specific heat measurements on three samples of polyethylene, 
differing only in density, were made in the temperature range of 2.4 
to 30°K. A definite density dependence was noted for the specific 
heat in this temperature interval which allowed extrapolation of the 
data to completely crystalline and completely amorphous cases. At 
the lowest temperatures, the amorphous results were observed to dis- 
play an "excess" heat capacity which could be accounted for by the 
occurance of a single delta function peak in the low frequency part 
of the vibrational spectrum. This excess did not appear in the com- 
pletely crystalline extrapolated data, and the specific heat was found 
to be proportional to the cube of the temperature up to 99K. 

An attempt was made to compare the results with two theoretical 
models (Tarasov and Stockmayer-Hecht) with only fair agreement in one 
case, and none in the other. The agreement with previous experimental 
results is excellent if the density dependence and the excess heat 


capacity are considered. 
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1. INTRODUCTION 
General Remarks 

The objective of the investigation which is the subject of this 
thesis was to determine the specific heat of polyethylene by direct 
calorimetric methods in the temperature range between 2.4 and 300K. 
This study was prompted by some earlier results of low precision which 
indicated a rather sizeable density dependence of the specific heat in 
the low temperature range [34] . This observation was in contradiction 
to statements in the literature which denied such a dependence [52] . 
In addition, theoretical models had predicted an unusual behavior of 
the specific heat of fibrous solids, such as polyethylene, at low 
temperatures which had not been verified by previous measurements 
[19, 46, 47, 52 . it was hoped that these effects, if they were to 
occur, would be found in the range covered by the present experiment 
thus affording a comparison of the experimental results with the pre- 
dictions. Further, it was desired to investigate the density depend- 
ence over a larger temperature range, and with more precision than 
that in the work noted above. 

The unusual heat capacity behavior which had been predicted to 
take place is a consequence of the polyethylene molecules being long 
chains, and the forces along the chains being much larger than those 
between chains. This one dimensional characteristic gives rise to 
theoretical descriptions for the specific heat which differs from 
most other non-fibrous solids. The theoretical models based on this 
chain structure will be discussed later, as will the morphology of 
polyethylene. However, for now it should be noted that polyethylene 


is an example of a polymer which crystallizes partially, and thus is 
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neither completely crystalline nor completely amorphous. So, one 
parameter in any discussion of polyethylene must be of its crystal- 
linity, which is directly related to its density. Making the assump- 
tion that polyethylene consists of two distinct forms, the crystalline 
and the amorphous, it is much simpler, for this discussion, to des- 
cribe the specific heat of polyethylene in terms of the specific heat 
of these two idealized phases, i.e., the specific heat of any real 
sample can be expressed as an admixture of that of the two hypotheti- 
cal ones. In order to accomplish this, the heat capacity of various 
samples differing in crystallinity (density) was determined so that 
extrapolation to the completely amorphous and crystalline cases could 
be made. 

Heat capacity data for polyethylene is available from various 
sources in the temperature ranges from 1 to 59K and from 209K to its 
melting point (about 428°K). The results of this research will 
bridge the existing temperature gap and should be of interest for two 
reasons. First, since polyethylene is an archtype of the chain 
polymers, the behavior of its heat capacity over the entire tempera- 
ture might be representative of polymers in general, and second, to 
allow a possible check on the aforementioned theoretical predictions. 

Previous Measurements 

Below 100°K there is not much experimental data on the specific 
heat of polyethylene. Above this temperature a considerable number 
of measurements have been made to determine the specific heat of the 
solid as it passes through the glass transition and the melting point. 
However these latter results are not germane to the present investiga- 


tion. 
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The data for the temperature range below 5°K were available from 
two sources. Isaacs and Garland reported measurements on a single 
high density sample in the 1.8-5.3°K interval obtained by calorimetric 
methods moms Ihe heat capacity of three samples of various densities 
was measured by indirect means in the range of from 1 to 4.5°K by 
Reese and Tucker [34]. The density dependence noted before is quite 
apparent when all four sets of measurements are compared [34]. It 
should be noted, however, that the data was analyzed assuming a cubic 
temperature dependence for the specific heat. This assumption is not 
verified by the present investigation, which shows a distinct non- 
cubic relation in this temperature region; however, the density de- 
pendence is valid. 

Dainton, Hoare, Evans, and Melia give results for two samples in 
the range 20-3009K [13] . However no conclusions could be made about 
the density dependence, since the samples were very similar in density. 
Ihe only other known source of specific heat data for polyethylene in 
this temperature range is the work of Sochava from 17 to 60°K [un]. 
The sample used by Sochava was not well characterized, the only infor- 
mation given being that it was a "typical commercial product", Thus 
there is no way to determine the density relatio@ships in the higher 
temperature range above 20°K. Above 90°K there are more measure- 
ments but it is felt that this information was not of interest insofar 
as the present investigation is concerned. 

Brief Description of Polyethylene 

Chemically, polyethylene is one of the simplest forms of chain 

polymers, consisting only of CH, units. The bonding along the chain 


is covalent. The chain length can be varied during manufacture, but 


1. 


in general it is quite long, on the order of 1000 units. In addition, 
side branches can and do form on the chain, and is also a parameter 
which can be varied according to the manufacturing process. These 
branches, which are short chain sections normally only 2-5 units 

long, join to the main chain by the substitution of a carbon-carbon 
bond for a carbon-hydrogen bond [31,38] . 

Polyethylene is one of the polymers which crystallizes partially 
in the bulk form. A necessary condition for crystallization is geo- 
metrical regularity of the molecular structure, which is fulfilled by 
the simple form of the polyethylene chain Lie, 31]. This regularity 
allows the chains to form ordered arrays, uia are considered the 
crystalline state. However on the basis of single crystals grown 
from solutions, which are quite small, it appears that the crystals 
are lamella, approximately 100 Angstrom units thick, with the molecu- 
lar chain axis being perpendicular to the face [23,35] . Whether this 
crystalline form is the same in the bulk is not known conclusively [31] ° 
However, it is well established that the crystalline part of poly- 
ethylene is composed of clusters of spherulites [31] š 

The most important factor inhibiting crystal formation is branch- 
ing. This is evidenced by the fact that highly branched polyethylene 
displays relatively low crystalline content [45] . This is quite 
reasonable since the regularity of the structure is broken at the 
branch points, thereby producing discontinuities or distortions in 
the crystal. 

If the chain is not in the well ordered array of the crystalline 
state, it is in some metastable position, intertwined and twisted 


around other chains. The complexity of the non-crystalline form makes 
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description very difficult. However, it sufficient for this discus- 
sion to characterize the non-crystalline part of the polyethylene as 
the amorphous regions. 

Since in the crystalline part of the polyethylene the molecules 
are more closely packed than in the amorphous, a higher density is 
found for the crystalline phase than for the amorphous [15] — TI 
fact if the density, e ; of a particular sample is known, its weight 
percent crystallinity, ^ » can be calculated from the relation, 

x. fh < 

1.1 ٩ Cc = CA ) 
where (c and (^ refer to the densities of the crystalline and 
amorphous parts which are abtained by indirect methods [15] . The 
erystallinities of bulk polyethylene vary in general between 35 and 
95%. 

Heat Capacity Models 

The specific heat of a solid can be calculated if the tempera- 

ture dependence of its internal energy,U, is known. This follows from 


the definition 


QU 
2 G. J Tr 


If the internal energy's temperature dependence is due to vibrations 
of the constituent parts of the solid about their equilibrium posi- 
tions, the total energy of oscillation over the whole range of the 


vibrational spectrum is, 


1.3 


A = (yo Ev, 7)4r , 
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Av/kT _ 
where رت‎ TS Aa Ke Ë is the average energy of a Planck 


oscillator with vibrational frequency Y, and 4 (r)is the distribution 
Function for the eigenfrequencies of elastic vibrations of the solid, 
which is called the frequency spectrum [ 14] . If there are N mass 
points in the solid, then there can be only 3N normal modes of vibra- 
tion, assuming that the temperatures are sufficiently low that possible 
internal motion of the mass units can be neglected (the optical modes 


of vibration). This sets the restriction 


1e 5N = (gody : 


Rigorous treatment shows that the upper limits in the integrals in 
equations 1.3 and 1.4 can be replaced by a maximum or cut-off fre- 
quency Yi, [ 50 ik Hence if 4 Qr) is known, X. can be determined from 
equation 1.4. 

This cut-off frequency can be interpreted as the frequency at 
which two neighboring mass units differ in phase by 180°. It should 
also be noted that Vy, is related to the forces acting between the 
neighbors, that is, greater forces imply larger cut-off frequencies, 
assuming particles of equal mass [14]. 

Ihe central problem in the calculation of the internal energy 
of a system, and therefore its specific heat, is to find the frequency 
spectrum eo). This is in general extremely difficult. Therefore, 
several approximation methods have been employed for solids. One of 
the most successful approximations is the Debye model. 

The Debye approximations are that the frequency spectrum is that 
calculated assuming the solid is an isotropic continuum, subject only 
to the condition stated above concerning the cut-off frequency, and 


that there is a linear relation between the frequency and the wave 
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number, i.e., that there is no dispersion. These approximations can 
be expected to hold for long wavelengths (low temperature) [ 50] : 


Debye theory than predicts a frequency distribution given by 


9M ——* vey کله سر‎ 


W اھ‎ A miu. 


Using this in equations 1.2, 1.3 and 1.4 leads to the familiar Debye 


equation for the specific heat, 


3 e Us (#2) 
TN Op /T 


I2 CE- =) 7 RET — 3 CHEIT = | 


١ 


1.6 Ge 


۱١ 


where D; (a) 


y^ e ° ° e 
and $,- ga is a characteristic temperature of the material called 
se L 
the Debye temperature. For very low temperatures, 22 << 30? 


equation 1.6 reduces to, 





3 
ATi EE 
1.7 Cy= 28 -F T Tano 


and for very high temperatures, 


1.8 Cy, = $K T > ceo 


This last result is the same as is predicted classically using the 
equipartition of energy [24]. 

As a logical extention of this model, it is more reasonable to 
assume that there exists a cutoff wavelength, rather than a maximum 
frequency. Since it is known that for all solids the velocity of 
sound is dependent of the direction of propagation, the cut-off fre- 
quency will also be different for the different modes and directions 


[ 43]. Recognizing this, equations 1.5 and 1.6 become, 
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4 
کے‎ 
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x 
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1.9 oy (y)dy = LN Y 


H 2 Ve < YA, 
$ vs ¿ š 
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where the subscripts refer to the transverse and longitudinal modes, 
and the Y's are directly related to the average sound velocities in 
these two modes. It should be noted that this modification preserves 
the T? dependence of the specific heat at the lowest temperatures. 
Polymers differ from most other solids in the sense that they 
are essentially one dimensional in nature, because the forces actis 
along the chains are much greater than those between chains. If a 
Debye-type analysis is made for a one dimensional system, in a manner 
completely analogous to the preceeding, one obtains the following re- 


sults [ 47); 


1.11 o = ر‎ 4 y 
C) ۸ 
oi 
1212 C, 2 3R D, (S) 
A Ba T. d e 
where p (Ee 2 x) | — : ۶۹۸۲م‎ -١ 


4 


and O = Á Vo 


The model employed above describes non-interacting chains, so 
is not of direct interest in this discussion. However, if the chains 


are intertwined and linked with one another, as in a chain polymer 
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Such aS polyethylene, there will be both the three dimensional and one 
dimensional continua existing together. This is the basis of a model 
due to Tarasov [ 47 ] . Since the forces acting along the chain are much 
greater than those between chains, it is expected that the maximum 
frequency for the one dimensional part of the model will be greater 
than that for the three dimensional part. Tarasov therefore predicted 


that there would exist two frequency spectra such that, 


7 دلا 
3N = f d (w)dw + & (v) dy‏ 1.13 
i 1‏ 


where the subscripts refer to the one and three dimensional quantities. 
Using the Debye assumptions, frequency distributions can be found. 
However, the separation into transverse and longitudinal parts is 

again applicable, as it was above. For the remainder of this deriva- 
tion only the longitudinal part will be considered, and the transverse 
contributions put in at the end by analogy with the Debye case.  There- 


fore equation 1.13 becomes, 


£ 4 
1.14 N = f A3) dv +f مین یئ‎ 
0 Jose 


for the longitudinal mode only (denoted by the superscript 2 ye 59 


equation can be rewritten as 


1.15 N = Nae FMS 


where N30 and N 3 refer to the first and second integrals above. 


Ihe distribution of the modes are, for N13 


/ dv 4 £ 
1.16 y rdv = TENIS y, 4- y; Y>; < y < ,لا‎ 
O VOV y < yz 
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and, for Nao, 


۾ y-‏ 2 پس v‏ ہے 
qi (v)dw - A v Os) d Í‏ 02 
٥ y. I‏ 


If these spectra are the correct ones, the heat capacity is given by 


the resulting Tarasov equation, 


ps J‏ - )3( 0( - - زک 


1.15 


for the longitudinal mode, where Gi Š, are related to وو‎ 
and the other symbols have been previously defined. Inclusion of 


the transverse modes then yields the complete Tarasov relation, 
p L 
1.19 Cue Cy tace. 


In the two temperature extremes, equation 1.19 reduces to 


۱ 
nk Nn E =) 4 
1720 C, 23R E F کې‎ 934 - EIS O 


3K س‎ oO, 


ft 
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Between these extremes, there will be a region for which, 


1 2 l 

12 c,» se T [S Iz * m) T, 8 66%. 

No continuum model should be expected to explain the specific 
heat in the presence of dispersion, which will take place in any real 
lattice. In order to calculate the density of states in the vibra- 
tional spectrum correctly, the effects of all of the forces acting 
between the mass units must be included. In practice this is not 
feasible. As a result, the usual method of approach is to make real- 


istic estimates of the strength of the forces, and drop all but the 
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dominate ones from the calculations. Blackman, and many others since, 
have achieved a certain amount of success along these lines; the de- 
gree of success depending on the validity of the assumptions about 

the relative strength of the forces, and the number of forces included 
in the calculations [ 6,7,16,28 ]. This type of analysis, since it em- 
ploys the discrete nature of the solid, should be more complete than 
the continuum treatment, if the correct parameters are used. There is 
however one feature of the lattice dynamics calculations which is the 
same as that obtained by the continuum assumptions: the heat capacity 
at the lowest temperatures should be proportional to the cube of the 
absolute temperature. 

For crystalline chain polymers an attempt has been made to obtain 
the vibrational spectrum using lattice dynamics by Stockmayer and 
Hecht whose work was greatly extended by Genensky and Newell [19,46]. 
This analysis is based on the consideration that there are strong 
forces along the chains and weaker ones between them. The frequency 
spectrum deduced was quite complicated, but has certain characteris- 
tics which are predicted to be characteristic of the heat capacity of 
fibrous solids in the crystalline form, and thus should be applicable 
to crystalline polyethylene. These predictions include, (where Cj 
now refers to the specific heat contribution from the interchain 


interaction, and C3 to that along the chain) [19] 





4C,4 C4 
BS Cy = — او‎ 
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«T, 
1.28 Ca په‎ 
where [^ is proportional to the maximum vibrational frequency VQ, 
c is related to the nearest neighbor forces, and K to the forces which 


oppose the bending of the chain. The approximations are, 
1.29 eras B/. 2 ٥ ) (o9 


where P is related to the force constant along the chain, and ۷ is 
related to the second nearest neighbor forces. Using these relations 
temperaturei] regions are predicted for which the heat capacity is at 
first (lowest temperatures) proportional to the cube of the tempera- 
ture, then changes to a TP dependence, where b lies between 2.5 and 
3, than b decreases to approximately .8 for T about one half s 

To complete this brief treatment of specific heat models which 
might be applicable to polyethylene, a related phenomenon, which 
appears in the low temperature heat capacity of some amorphous solids, 
must be mentioned. If the heat capacity of these substances is 
measured calorimetrically at low temperatures, an "excess" specific 
heat is observed. To understand the term excess one should refer to 
equation 1.7 which can be rewritten as [ 24], 

تر ے om‏ 


a ee 


1.30 v 7 (en FE ir 3 


where V is the volume of the sample, and ùv is the velocity of sound 


at these temperatures, which can be interpeted as, 
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l | = x ET < 
لے شارت ےت وق‎ 
Therefore sound velocity measurements should yield the same value for 
C/T? as those obtained from direct calorimetric methods. However, it 
has been found experimentally that for some amorphous substances, such 
as vitreous germanium, vitreous silica, and glycerol giass, the calor- 
imetric results are much greater than those given by the acoustic 
methods and are not T2 dependent Ip» 17 |. 

The explanation for this anomolous effect is that at low temper- 
atures there may be modes excited which are not measured by the 
acoustic methods. Since the mechanical sound measurements are con- 
cerned with the averaged (in the sense that the measurements must be 
macroscopic) elastic forces between the particles making up the solid, 
if there are low frequency modes which do not enter into these elastic 
vibrations, then the calorimetric measurements should be higher than 
predicted by equation 1.30. To explain the existence of these extra 
modes, Rosenstock has suggested that in disordered lattices, such as 
the amorphous solids, microscopic holes exist which may contain a 
particle which is bound in only one or two directions [ 39,40 ]. These 
loosely bound mass units can be likened to particles free to vibrate 
within a definite volume (the hole), and could thereby account for the 
low frequency non-acoustic modes. However, there are other hypotheses 
which attempt to explain the modes as being due to other phenomena 
such as non-symmetric lattice sites or distorted bonds [ 7k In 
any case, the origin of these low frequency modes is not clear at 
this time. In the case of the examples cited above it has been possi- 


ble to empirically describe the low temperature departure of the 
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specific heat from the T? relation by adding a term (or terms) pro- 
portional to the Einstein specific heat to equation 1.30 آ0‎ 
This implicitly implys that some of the vibrating units do not con- 
tribute to the acoustic odes and are accounted for in the total 


specific heat by an expression such as, 


2 Cex - 3 = A E f) 


وم 


Aw Se) = e ue ) pu 


where Or = E 967ب‎ = 


Ve is the excitation frequency for these modes, and Mz is the num- 
ber of vibrators involved [24]. Adding a term such as this to equation 


1.7 or 1.30 results in 


e C, = AT +E E (FE) 


where A and [3 are constants. Since the specific heat is no longer 
proportional to the sube of the temperature, a plot of the C/T? versus 
temperature will demonstrate the presence of the excess heat capacity. 
2. EXPERIMENTAL DETAILS 

Introduction 

In all calorimetric investigations it is necessary to know two 
quantities; the temperature of the sample and the acount of heat added 
to it. In this experiment, meeting the former requirement was the 
more difficult. 

The temperature determination must be made quickly and accurately. 
Since resistance thermometry fulfills these requirements, doped ger- 
manium resistors were selected as thermometers. These were chosen 
rather than carbon resistors since germanium resistors are reputed 


to maintain their resistance-temperature relationship even after a 
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large number of cyclings between room temperatures and liquid helium 
temperature [18,25] . Other experimenters have found that this is not 
true for carbon resistors, so that they must be recalibrated after 
each temperature cycling [18] . The need for recalibration would 
greatly complicate the experiment thus making the choice of germanium 
resistors an easy one. 

Since the temperature range covered by this investigation was 
from about 2 to 309K, doped germanium resistors (Minneapolis-Honey- 
well, models 2401 and 3401) were selected which had sufficiently high 
sensitivity, and yet, for the above temperatures, had resistance 
values within the range of ordinary laboratory resistance measuring 
devices. Once the resistors were chosen the most important task was 
to find the relation between the resistance and the absolute tempera- 
ture. The discussion of the calibration procedures is the subject of 
Appendix A. 

Having insured that a quick and accurate determination of the 
sample temperature and its changes could be made, attention was di- 
rected to the second requirement. Adding heat electrically to a 
resistance element in thermal contact with the sample provides a 
simple and accurate method of deliberately raising the sample's tem- 
perature. However great care must be taken to insuré that this 
energy goes solely to the sample, and that it is the only source of 
heat input. To this end, an adiabatic calorimeter was constructed in 
which the sample could be suspended in vaccuo by high thermal resis- 
tance mountings from a heat shield, which completely surrounded the 
sample and was held at the same temperature. The design and construc- 


tion of this calorimeter will be discussed in a subsequent section. 
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Sample Preparation 

The heat capacity of three polyethylene samples of various den- 
sities were to be measured. Two of the samples were Marlex, manu- 
factured by Phillips Petrolium Company, and were received in the form 
of extruded cylindrical pellets approximately 3 mm. in diameter and 
length. One of the samples, Marlex II, had branching introduced to 
impede crystallization. The third sample, a typical commercial low 
density polyethylene, was a 1" round bar. The samples are referred 
to as Marlex I, Marlex II, and Low Density in the remainder of this 
thesis. 

An important consideration in the preparation of the samples was 
that the time interval between the application of heat and the achieve- 
ment of thermal equilibrium be as short as possible so as to minimize 
errors due to heat leaks into or out of the sample. The low thermal 
conductivity of polyethylene at low temperatures. makes this a diffi- 
cult problem [ 34]. Because of the calorimeter design, the samples 
were to be cast into 9/16" rods about 6 cm. in length. If such a 
polyethylene rod were heated from the center, it was estimated that 
it would require times of about 100 seconds to achieve thermal equili- 
brium at 5°K and nearly 1 1/2 hours at 30°K. Such comparatively long 
thermal relaxation times would undoubtedly introduce serious errors 
in the determination of the amount of heat added to the sample unless 
exceptional precautions in the control of the temperature of the adia- 
batic shield were taken, in addition to making the experiment extreme- 
ly tedious to perform. Thus san eee effort was expended to pre- 
pare samples with shorter thermal equilibrium times. The method 


chosen was in inbed within the sample 99.999% pure copper wires such 
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that the wires were no farther than .3 cm. apart. In addition the 
sample was wrapped in a copper foil. The heat could thus be distri- 
buted throughout the sample by means of low thermal resistance copper 
paths. This was expected to greatly reduce the thermal relaxation 
times estimated above. The success of this method is seen by the fact 
that the time required during the actual experiment to achieve thermal 
equilibrium was never more than about 12 minutes at the highest temp- 
eratures. 

Introduction of the wires into the sample in a uniform fashion 
presented a considerable problem. The following system was eventually 
evolved. The pellets (or in the case of the low density sample, the 
shavings obtained from turning the bar on a lathe) were first cast 
into small rods of 1/4" and 3/16", and wrapped with 30" of the copper 
wire. The rods-were compactly assembled within a 9/16" mold, and 
additional wires inserted into the resulting gaps. The sample was 
then cast at 155°C under 900 PSI for 16 hours, and then allowed to 
cool to room temperature (about 12 hours). 

The extra material added to the samples was to be considered 
part of the calorimeter's sample holder. It consisted of 30' of 
.0051" diameter drawn copper wire, .001" thick copper foil 5.580 by 
5.035 cm., and .08 grams of G.E. 7031 varnish. 

Ihe polyethylene density was determined by measuring the density 
of the entire sample and correcting for the copper wire which was 
introduced. The crystallinity of each sample was then calculated by 


use of equation 1.1 with the values 
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where the (^ and (e refer to pure amorphous and crystalline poly- 
ethylene densities [ 9,36 ]. It must be noted here that the densities 
given above are not the only ones available. For the amorphous case 
there is fair agreement in the literature that .85 is a good average, 
but more attention must be given to the crystalline density [ 15]. For 
this case, the density has been measured by X-ray methods to be as 
high as 1.014 gm/cm? and as low as .964 gm/cm? depending on the amount 
of branching present [ 49 ]. Furthermore voids or free volume in the 
spherulitic crystalline form of bulk polyethylene could reduce the 
higher figure noted above even in the absence of branching [ 10,11]. 
These considerations will be referred to again when the error of the 
experiment is discussed. The measured densities and calculated crys- 
tallinities are given in the table below, where the error indicated 
in the latter quantity is that due to the uncertainty in the density 
measurement only. 

TABLE I 


Densities and crystallinities of the samples. 


Density % Crystallinity 
Marlex I .973 = .002 gm/em3 ۰ 
Marlex II ء‎ 2 Tete Sie Tee WG 
Low Density -915 2 45.9 * 2.5 


Calorimeter Design and Construction 
In order to accurately measure heat capacities, the sample must 
be thermally isolated to insure that the noted temperature differences 
of the sample are due entirely to the measured heat input. To accom- 
plish this, it is necessary to remove all conduction and radiation 


paths to and from the sample. The method employed here was to enclose 
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the sample in a shield which is kept at the same temperature as the 
sample, but insulated from it. 

In this experiment the heat capacity would be measured at various 
temperatures above the temperature of a liquid helium bath by success- 
ive applications of heat to the sample. Hence the shield could not be 
in good thermal contact with the bath. This thermal isolation was 
accomplished by suspending the shield inside of a vacuum container 
which was in turn in contact with the bath. A can type adiabatic 
shield was made from a copper 2" tube, 3" long, threaded at the top 
to mate with its lid. The lid was soldered to a 3/8" stainless tube 
(.01" wall thickness), which connected to a suitable vacuum container. 
The tube then extended out of the cryostat to peripheral vacuum equip- 
ment and served as the pumping line and electrical lead housing. The 
shield also had provisions for receiving a resistor which would sense 
the shield's temperature. 

There would be a heat leak from the shield to the vacuum con- 
tainer through the support, but this was compensated for by means of 
an electrical heater on the shield. In addition this heater was used 
to raise the shield's temperature to that of the sample. Seventeen 
feet of .002' manganin wire, coiled and glued with G.E. 7031 varnish 
to the top of the shield, served as the shield heater. 

The sample holder, in addition to holding the 9/16" sample, was 
used to mount the sample heater, temperature sensing resistor, and a 
receptacle which would be the contact point for the thermal switch 
(to be discussed later). A copper ring with the necessary provisions 
was constructed for this purpose. When heat capacity measurements 


were made, the measured values would be the total of the sample plus 
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the addenda (the holder, heater, thermometer, etc.). Although the 
latter's heat capacity would be determined experimentally and sub- 
tracted from the total, it was still desirable to keep it as small 

as possible. Based on the size of the calorimeter, it was determined 
that a sample of 10 cm? could be accommodated. With a sample of this 
size it was estimated (correctly, as was determined later) that the 
heat capacity of the addenda would amount to less than 10% of the 
total, which was deemed acceptable. 

The sample heater, which consisted of 15" of .002" manganin wire 
(100 ohms), was would around the thermometer receptacle. A four wire 
connection was made to the heater, and the wires were thermally 
grounded to the adiabatic shield. The power was measured potentio- 
metrically using a Leeds and Northrup K~3 potentiometer. A Varian 
model 501A time mark generator, used in conjunction with a Hewlett- 
Packard model 521C electronic counter, determined the heating times 
(the timer was later changed to an Anadex model CF-200R counter- 
timer) and was triggered by the sample heater switch. 

It would be impossible, regardless of the care taken, to insure 
that the shield would be at exactly the same temperature as the sample 
at all times. However, if the shield were thermally isolated from 
the sample, small temperature differences could be tolerated. Based 
on this the sample holder was suspended within the shield by 6 nylon 
threads, approximately .12 mm. in diameter, and a minimum of 1.9 cm. 
in length. Assuming a maximum of .59 temperature difference between 
the shield and the sample, the heat leak through these supports would 
be on the order of 107? mj/sec which, for equilibrium times of less 
than an hour, could be neglected (radiation losses are of the same 


order of magnitude). 
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When the calorimeter was initially cooled, or when cooling back 
to the bath temperature after a series of measurements, a method had 
to be incorporated in the design to cool the isolated sample in a 
reasonable length of time. For this reason a mechanical switching 
device was constructed. This was chosen over the use of exchange gas 
in the vacuum spaces, since the latter method has been shown to have 
quite long recovery times and has attendent gas absorption problems 
Poo) The switch consisted of a brass finger which would make contact 
with a receptacle on the sample holder. The finger was soldered to a 
brass rod which extended out of the vacuum container into the liquid 
helium bath through a nickel bellows and from there out of the cryo- 
stat. The finger could be raised and lowered manually so as to make 
thermal and mechanical contact with the holder whenever it was re- 
quired. The original design was modified after an initial trial by 
tinning the brass finger and its mating receptacle with metallic in- 
dium, for better thermal contact. In addition a brass leaf was added 
which would make contact with the adiabatic shield when the thermal 
switch was closed in order to cool the entire calorimeter simultane- 
ously. 

The two thermometers (resistors) were attached to the shield and 
the sample. Three wire connections, which were thermally anchored to 
the shield, were made to each resistor. The resistance measuring 
equipment was identical to that used by Reese and Tucker and will not 
be discussed here [ 34 ]. 


The completed calorimeter is shown in figure ٠ 
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Figure 1, Calorimeter with samole in place, 
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Measurement Procedures 

In preparation for the heat capacity measurements, the sample 
was placed in the calorimeter and cooled to the desired dogs tem- 
perature with the thermal switch closed. When the thermometers indi- 
cated that the sample and shield were at the liquid helium bath temp- 
erature, the thermal switch was then opened. This set up a small 
amount of frictional heating as the contact was broken, but within 
30 minutes the system would again be in equilibrium. 

The sample was then heated by applying electrical power to the 
sample heater in amounts which would raise the temperature of the 
sample in 576 increments of its average temperature. The power was 
recorded at each temperature step as well as the initial and final 
readings of the time and the thermometers. The heating times were 
kept between 20 and 60 seconds, which required that the power be in- 
creased every four to six measurements. When the desired final temp- 
erature was reached, the thermal switch was again closed, the sample 
cooled, and the procedure repeated. 

uode the heat was being supplied to the sample, the adiabatic 
shield was also heated by means of its heater. This required care 
since, as previously mentioned, it was highly desirable that the 

, 
shield and the sample be at the same temperature. After several 
trials, however, certain regularities in the amount of heat needed 
for the shield at various temperatures became evident, and the desired 
result was obtained. It is felt that there never existed more than a 
one percent temperature difference between the two thermometers while 
actual measurements were being taken, and for the bulk of the runs, 


the difference was considerably less. 
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Due to the close proximity of sample thermometer to the heater, 
while the power was being supplied the temperature would be observed 
to overshoot the desired value, and, when the power was removed, drop 
back exponentially. The time required for the sample to reach equil- 
ibrium was from 2 seconds to about 12 minutes as the temperature was 
increased from its lowest value to its highest. These times were not 
considered excessive, except at the highest temperatures when the 
shield might tend to drift while equilibrium was being established. 
This drifting was countered by monitoring both thermometers and keep- 
ing the shield at the desired temperature by removing or increasing 
the shield power. 

The heat capacity of the samples was determined by five separate 
sets of measurements. ‘The first three covered the range of 4.5 to 
30°K, one series for each of the three samples. The other two were 
for the temperature interval 2.4 to about 89K for the Marlex I and 
Low Density samples only. For each set of measurements the aforemen- 
tioned temperature ranges were covered three times. 

Since the measured heat capacities were the sum of those of the 
sanple and the addenda, the latter was measured in two separate ex- 
periments, for the temperature ranges of 4.5 to 300K and 2.4 to 89K, 
so that its contribution could be subtracted from the total. In each 
of the experiments, the addenda heat capacity was measured four times 
over the covered range. The procedures were identical to those des- 
cribed above except due to the reduced mass, the amount of heat re- 
quired and the equilibrium times were both substantially decreased. 
3. RESULTS OF THE EXPERIMENT 

Data Reduction 


When the experimental measurements of the heat capacity were 
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made, the data consisted of the power reading, the elapsed time inter- 
val during which the power was applied, and the initial and final 
temperatures as indicated by resistance values. The heat capacity 
can then be calculated from, 

3.1 C = — Power X Elapsed Time , 

TPinal -- 2۱1 

Due to the complexity of the resistance-temperature relation (Appen- 
dix A), the calculations were performed by the CDC 1604 digital com- 
puter. 

All of the data were reduced to heat capacities, and of the ones 
which displayed large deviations, only those attributable to obvious 
blunders in data recording were rejected. 

Since the heat capacities measurements for each sample represen- 
ted the total of the sample plus the addenda, the latter's contribu- 
tion had to be taken into account and corrected for. As mentioned 
before, measurements of the heat capacity of the addenda were per- 
formed. These heat capacities were for discrete, random temperatures 
and, as such, were not useful for the adjustment which was to be made. 
Therefore the data were best fitted, in the least squared sense, by 
the computer to a continuous heat capacity-temperature relation of 


the form, 
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With this relationship the total heat capacity for each sample 
could be corrected so that the specific heat of the sample alone 
could be determined. These calculations were also performed by the 


computer and the resulting heat capacity data were again fitted to 
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equation 3.2. The specific heat was obtained by dividing the heat 
capacity by the mass of the sample. 

The measured specific heat for each sample is found in Appendix 
C, and the parameters of equation 3.2 for the samples and the addenda 
in Appendix E. The smoothed results from the fitted curve are shown 
in Appendix D, and a specific heat-temperature graph is displayed in 
figure 2 where the difference between the smoothed curve and the actual 
data points is too small to be shown. 

The temperature range, total number of data points, and the per 
cent RMS deviation of the actual data from the smooth curve are pre- 
sented in Table II. 

Table II 


Temperature range, amount of data, and RMS fit of the data 
to a smoothed curve for each of the samples and the addenda. 


Temperature Number of HRMS 
range data points deviation 
Addenda 2.4 - 300K 218 1.70% 
Marlex I 2.4 - 30 161 1:552 
Marlex IT 4.5 - 30 Iob 1.01 
Low Density 2.4 - 0 158 1.09 


Extrapolation to 100% and 0% Crystallinity 
As mentioned in the first chapter, the specific heat of any 
polyethylene sample might best be expressed as admixture of the specif- 
ic heats of the completely amorphous and completely crystalline phases. 
As suggested by the earlier low precision work of Reese and Tucker, aad 
very well substantiated by the present investigation, the relation is 


linear [34]. That is 
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where the C's on the right side refer to hypothetical values for the 
specific heat of crystalline and amorphous polyethylene. This lin- 
earity is reasonable if one considers that it is possible (in theory) 
to break a sample up into many very small macroscopic volumes, each 

of which is either entirely amorphous or crystalline (if the boundar- 
ies are neglected). The two idealized polyethylenes should have dif- 
ferent values for the specific heat at any particular low temperature 
since the amorphous state should be characterized by a set of para- 
meters which will differ from those of the crystalline state inasmuch 
as the interchain binding scheme is completely different. 

Using this assumption, the heat capacities for the three sam- 
ples, divided by the cube of the temperature, were plotted as a 
function of crystallinity for various temperatures (Figure 3). The 
linearity is quite apparent, and on the basis of this, it is felt that 
equation 3.3 is indeed valid. Using this linear relationship one can 
perform a simple extrapolation to obtain values for the specific heat 
as a function of temperature for the hypothetical 100% and 0% cases. 
However, rather than use the C/T relation, which is shown for des- 
criptive purposes only, the computer was used to perform the extrapo- 


lation from the equation 


3.4 C s X 4 


where /Y^ and Ar were determined from the experimental (smoothed) 
values of C, the specific heat for the three samples of differing 
crystallinities. The extrapolation to both 0% and 100% was performed 
every .049 and the results were fitted to equation 3.2, the same as 
the real samples. 


Appendix D presents the extrapolated data, and figure 2 is a 
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Figure š, C/T) versus crystallinity. 
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plot of the specific heat versus temperature. In Appendix E are 
shown the parameters for equation 3.2 for all of the real and extra- 
polated samples. 
The 100% Extrapolated Data 

The heat capacity obtained from the extrapolation to 100% 
crystallinity had some interesting features. A plot of C/T versus 
temperature shows a relatively large temperature region (up to about 
99K) where the specific heat has a cubic temperature dependence (fig- 
ure 4). This is in agreement with the expectations if there exist 
relatively large binding forces in the crystalline case, since on the 
basis of the Debye model larger forces imply larger @ § which in turn 
imply an increase in the temperature range in which the cubic depen- 
dence holds (equation 1.7). 

For this case, an excess heat capacity is not indicated by the 
C/T? plot, nor is one expected since this material should be a con- 
tinuous polycrystal. 

When the heat capacity for this hypothetical sample is plotted 
against temperature on log-log paper, it is possible to obtain the 


parameters A and a of the expression, 
o. 
3 5 CEAT, 


There appeared two regions in which they were constant, indicating 

that the specific heat is proportional to some power, a, of the temp- 
erature. The first region was the cubic dependence noted before. The 
other was from about 20 to 26°K where the exponent, a, was approximate- 
ly 2.37-2.40. In between the value decreased uniformly. This is in 
disagreement with the prediction of Genensky and Newell who expected 


it to be equal to some constant slightly greater than 2.5 for an 
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Figure 4, 
A plot of C/T% versus temperature for polyethylene, 
where the solid lines are the experimental and the 
extrapolated data, the dotted line below the 0% curve 
is the data minus the Einstein term, and the other 
dotted lines are the vredicted T z O°K values for 
these samples, The dots represent the predicted 


specific heat values in all cases. 
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appreciable temperature region, for a crystalline polymer [19 ]. 


This contradiction will be discussed later. 

The Tarasov model could be well fitted to the extrapolated data. 
The method used to determine the parameters was essentially the same 
as that suggested by Wunderlich, except that the equations used are 
for the extended model, equation 1.19 rather than the one suggested by 
farasov, which doesn't separate the longitudinal and transverse con- 
tributions [ 52]. According to equation 1.22 there should exist a 
region in which thé heat capacity becomes linear in temperature. This 
will occur at a temperature well above that at which most of the three 
dimensional modes are filled and above the temperature range of this 
investigation. Therefore the data of Dainton, et al, giving the 
specific heat of polyethylene from 20-3009K, were used to determine 
the linear region iu. This procedure will be justified later. The 
linear region appears to occur between 80 and 110°K, where C/T has 
the value 6.69mj/9K?-gm. Rewriting equation 1.22, it is seen that in 


this region, 


ps — + (| 
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Another relation between Cir and fie is necessary to determine them 
uniquely. This was accomplished by considering that the Oi, since it 
was the longitudinal mode limit for the chain, might be well character- 
ized by the largest frequency for which the C-C bond will be stretched. 
Reference to the Raman spectrum of polyethylene indicates a value of 
1150 cm”! for this parameter, and therefore a value of 1655°K was 
assumed for 1. [15]. This fixed Gir at 707°K. (The use of an appar- 
ent optical mode in this discussion, which up to now had only consid- 


ered acoustical vibrations, will be discussed later). 
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Having determined these two parameters, it is then necessary to 
fix the remaining two thetas needed for the extended Tarasov model, 
zr and aL. One relation comes from equation 1.20, which is here 


rewritten, 
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This is the low temperature limit, where C/T is constant, and the re- 
sults below ^~ 9°K could be used to evaluate the left hand side. In 
order to solve for either 4r or Os, another relation is needed. 
This was obtained from sound velocity measurements at liquid nitro- 
gen temperatures for a highly crystalline sample. These measurements 
were performed by Professor W. Reese [ 33 ]. Extrapolation of the re- 
sults to zero degrees Kelvin predicted that the ratio of the longi- 
tudinal sound velocity to that of the transverse was 1.67. Since the 
thetas have the same relative relationship to one another as do the 
sound velocities, Üsr and 3, could now be determined. The Tarasov 
parameters for this case are found in Table III. 
The O% Extrapolated Data 

Extrapolation to the 0% crystallinity (completely amorphous) 
case led to results which were quite different from those discussed 
above. In particular, the plot of C/T? versus the temperature shows 
that the heat capacity was not proportional to the cube of the temp- 
erature anywhere. The curve starts with a positive slope at 2.59K, 
increases to a maximum at around 59K, then decreases uniformly. This 
is also the case for the real samples as can be seen in figure 4. 
Since even the Marlex I sample (crystallinity of 84.1%) displays this 


hump, while the 100% extrapolated sample does not, it is quite apparent 
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that the presence of the hump must be due to the amorphous poly- 
ethylene. The hump is interpeted as an indication of an excess heat 
capacity. Referring to the first seetion the heat capacity for the 


lowest temperatures should then be of the form, 
= I 
200 ea Go 
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where í, A and Ü£ are constants to be determined. This equation 
was best fitted, in the least squared sense, to the extrapolated re- 
sults (and to the Low Density data where a correction was applied to 
account for the crystalline part) to yield values for Y, A and O¢ . 
For the model which describes the portion of the heat capacity not 
given by the Einstein term, one is forced to accept a continuum type, 
since it is not possible to do lattice dynamics for the amorphous 
case. From the foregoing discussion, the extended Tarasov model was 


chosen. Assuming this, A is interpeted as 
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There is no reason to assume that the one dimension parameters 
will change from those used previously, since they are merely func- 
tions of the forces along the chain. However another relation be- 
tween D and OL is still needed in order to determine them uniquely. 
The ratio of the longitudinal and transverse sound velocities used 
earlier is not really applicable in this case, since the sound vel- 
ocities were those for a highly crystalline sample; but, lacking any 
other evidence, the same relation was assumed. 


Using the above assumptions, and the values of A and X , the 
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Tarasov parameters can be found. They are presented in Table III 
along with othér results for this extrapolated hypothetical sample. 
Summary 

The percentage difference between the 0% and the 100% extrapola- 
ted samples were plotted against temperature as shown in figure 5. 
Although the curve is approaching zero, it is impossible to determine 
with any degree of certainty when if ever this would occur. However, 
it is felt that the percentage difference would be insignificant by 
1009K, where the data is linear in temperature [ 13 ]. This supports 
the assumption that €, will be the same for all samples. 

The heat capacities for the real samples and the extrapolated 
ones are shown in Appendix D at .59K increments, together with the 
specific heats predicted by the Tarasov equation, equation 1.19, which 


is here expanded as 
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The superscripts on the fs refer to either the 0% or the 100% para- 


meterg, which are presented in Table III. The crystallinities, X 5 
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are shown in Table I. 

Two hitherto undefined terms also appear in Table III. These 
are the Tarasov parameters which would result if no separation into 
the longitudinal and transverse modes were assumed, which is the ori- 


ginal model discussed by Tarasov. In this case, equation 3.11 would 


become, 
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In addition, the classic Debye temperature is shown as calculated 
from the low temperature limit. 

In figure 4 the dotted line beneath the 0% crystallinity sample 
curve is the result of subtracting the term Y (S5) cron the data, to 
illustrate how it accounts for the excess heat capacity. For the 
other samples, the dotted lines indicate the O°K intercept, and for 
all cases, the circled points are the values predicted from equation 
SO. 

In Appendix D the smoothed values of the specific heat for each 
sample is shown, along with the calculated specific heat of the Tara- 
sov model as given by equation 3.11. This calculation was performed 


by the computer (the mass of the polyethylene vibrating unit is taken 


as 14.1 atomic mass units). 
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Table III 


The Debye, Tarasov, and excess heat capacity parameters, 
and the low temperature values of the specific heat for 
the extrapolated data. 


100% 0% 

yj (m j/?K-gm) 0 . 3.24 
€ (°K) 260.1 175.4 
€; (°K) 873.9 873.9 
پټ‎ (OK) 141.9 78.6 
تی‎ KOK) 1655.0 1655.0 
مود‎ ٢( 707-0 707.0 
Oe, (°K) 222.8 123.4 
+ (OK) 227 75.1 
Y. (°K) - 28210 
C-Y£(5) aj/OKign, low .00787 .0256 
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Accuracy of the Results 

To estimate the accuracy of the heat capacity data, four sources 
of error were considered; the determination of the amount of heat in- 
troduced into the sample, the temperature determination, the effect 
of the shield, and the heat capacity of the addenda. 

The amount of heat energy introduced into the sample was obtained 
from three measurements; voltage, current, and time. The voltages (and 
current) could be determined to six places on the potentiometer. Since 
the same power was used on from four to six successive measurements, 
it was possible to check the reproducibility of the readings. Only 
at the very highest temperatures, where the resistance of the heater 
was probably changing with temperature, did it vary as much as .1%. 


The current measurement had the same inconsistency, but since it was 
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determined by measuring the voltage across a standard resistor, it 
also had an inherent maximum error of .05%. Timing was presumed ac- 
curate to one part in thirty thousand. The maximum total error from 
this source is therefore estimated to be .05% at the low end of the 
temperature s@@le and .25% at the high end. 

As discussed in Appendix A, the thermometers were assumed cor- 
rect to within .5% of the absolute temperature. This uncertainty will 
influence the accuracy of the heat capacity results in two ways. Since 
a temperature difference is divided into the total heat to obtain the 
heat capacity, an error can arise there. However any experimental or 
systematic errors in the temperature determination will be the same 
over a Short temperature interval and should cancel out (this is en- 
sured by the use of the smooth fitted curve which should average out 
these errors). Therefore the possible error from this source should 
be only that due to the error of judgment in reading the resistance 
bridge. This results in uncertainties in the heat capacity determi- 
0251711 E wile? 4°. upAbos3hout .279 at 309K; 

The other manner in which the inaccuracy of the temperature 
scale enters into the error is in assigning a temperature to a par- 
ticular heat capacity measurement. Since the specific heat can be 
expressed at any point as being proportional to some power of the 
temperature (the exponent will be between 3 and 2.4 for this case), 
the percentage error in the specific heat will be larger than that of 
the temperature. In this case the maximum possible error amounts to 
1.5% at the lowest temperature and 1.2% at the highest from this 
source. 


Since it was impossible to completely isolate the sample from 
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the shield thermally, considerable care had to be taken to ensure 
that they were at the same temperature, to prevent possible heat 
leakage. Experimentally, the effects of large temperature differ- 
ences were investigated, by purposely overshooting the shield tem- 
perature and observing the effect it had on heat capacity data. On 
the basis of these experiments it was decided that, with the experi- 
mental procedures used when actual data were being taken, this con- 
tribution to the error could be neglected for all but the highest 
temperatures, when, due to the insensitivity of the shield thermom- 
eter, difficulty insuring equal temperature of the shield and sample 
might occur. At these temperatures, the error might approach .5%. 

The heat capacity of the addenda was measured over the entire 
temperature range, so that it could be applied as a correction to 
the total measured for any sample. Since the mass of the addenda 
was chosen to be as small as possible, so that its heat capacity 
would be on the order of 10% of that of the samples, any possible 
errors which are present in the measurement of its heat capacity 
would be reduced by a factor of ten when applied to the total possi- 
ble error for a sample. Considering the preceeding this uncertainty 
amounts to approximately .2% over the entire temperature range. 

Considering all of the above, it is felt that the measured heat 
capacities were accurate to within 1.9% for the temperature up to 
about 26°K, and 3% above that. 

The subject of Appendix B is the measurement of the specific 
heat of copper. The results are compared with data found in the 
literature which presumably is of higher precision and accuracy than 


this work. The comparison shows that the results are consistant, 
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within the error estimate, but the present results tend to be low 
over most of the temperature range, which could be an indication of 
a ا تا‎ error. This point was not investigated further. 

As mentioned before, there was no statistical criteria estab- 
lished for rejecting any of the data, and the only points discarded 
were ones in which a large deviation could be attributed to a blunder 
in data accrual. 

When the data were fitted to equation 3.2, the RMS deviation 
(Table II) was found to be well within the estimated accuracy of the 
experiment and therefore the use of the smoothed data is considered 
to be a valid representation of the results. 

Since one of the objectives of this investigation was to verify 
the density (crystallinity) dependence of the specific heat of pé@ly- 
ethylene, and an extrapolation was made to the hypothetical 100% and 
0% crystalline phases, an estimate was made as to the reliability of 
this. The heat capacity of the three real samples would be exactly 
linear in crystallinity over the entire temperature range, if a varia- 
tion of at most 4% were made in any one of the heat capacities, or if 
a 1.5% adjustment were made in one of the crystallinities. Addition- 
ally there are other combinations of adjustments which could be made, 
but it is the uncertainty in the degree of crystallinity which is be- 
lieved to be subject to the most error. As noted in Table I, there 
is an inherent, approximate 2% error which is due solely to the un- 
certainties in the densities. In addition, since the densities of 
the amorphous and crystalline phases are not well established, it is 
estimated that the crystallinity of a sampl@ could be in error as 


much as 4%,ssbased on the density method of determining crystallinities 
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(equation 1.1). In view of this, it is felt that the linearity is 
valid to the precision of this experiment. 

Due to the possible error in the crystallinities, the accuracy 
of the extrapolated samples can be better than 2% for the amorphous 
case, and .8% for the crystalline phase. The uncertainties which are 
present in the actual specific heat is taken over into the extrapola- 
ted data. If the relative error in the heat capacity at any particular 
temperature is the same for all of the samples, as should be expected 
since they should exhibit the same systematic error if it is present, 
the maximum possible error in the extrapolated data due to this cause 
will be the same, or about 1.9% over most of the temperature range. 
Thus the uncertainty in these figures will be approximately 3 and 4% 
for the completely crystalline and amorphous phases. 

The precision of this experiment can be estimated by the close- 
ness of the fit of the actual data points to the, smoothed curve given 
- by equation 3.2. . The deviations appear random and the per cent RMS 


deviation is given in Table II, where the % RMS deviation is defined 


by, 





M t 
چا ء‎ 
3213 %RMS = 100% 
: N 
This precision, which is also an estimate of the reproducibility 
of the results, is between 1.01 and 1.52% for the polyethylene sam- 
ples. 
4. CONCLUSIONS 
Comparison with Other Results 
The results of this research have been compared with all of the 


known experimental work on the specific heat measurements on poly- 
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ethylene in their common temperature ranges. A detailed survey of 
these comparisons follows: 

(a) Isaacs and Garland [ 21] 

The specific heat measurements of Isaacs and Garland were on 
a single high density sample in the temperature range 1.8 to 5.3°K. 
They assumed that the specific heat in this range displayed T3 depen- 
dence and therefore did not separately determine the heat capacity of 
the addenda, but rather calculated it and subtracted it from the total. 
This led as they report to an unaccountable deivation from the T? be- 
havior, but they attributed it to experimental error. However it is 
felt that they were in reality observing the excess heat capacity 
which was disguised by their method of analysis. 

For the comparison in the mutual temperature range, it was noted 
that both the amorphous and crystalline parts of polyethylene are pro- 
portional to the cube of the temperature if the Einstein term is sub- 
tracted. Since their sample did not have the same degree of crystallin- 
ity of any of the ones measured in this work, the following equation 
was used to determine what specific heat should be expected on the 


basis of the present results; 


‘K 
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This is merely the low temperature limit of equation 3.12 with the 
measured values inserted. This assumed no model, however; only that 
the specific heat be T3 dependent at low temperatures, which it is 
for any realistic case, and that the linearity in crystallinity is 
valid. For the crystallinity of Isaacs and Garland's sample, two 


values were investigated; first the reported X-ray determined result 
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of 90%, and second, a crystallinity of 85% deduced from the reported 
density. In Table IV where the comparison is made, the two different 
values for the C/T3 resulting from the different crystallinities are 
subscripted to distinguish them. 

Examination of the table shows that the value deduced from this 
work with the density determined crystallinity is in excellent agree- 
ment with Isaacs and Garland's, if it is recognized that they could 
not separate the excess heat capacity. Use of this crystallinity 
rather than the X-ray value contains a certain amount of consistency 
since the present results were based on the density measurements. It 
is noted they give no estimate of the possible error in the X-ray 


value, although it is doubted that it is as high as 67. 


Table IV 


Comparison of present results with those of Isaacs and 
Garland. (Values expressed in m j/O K"- gm) 


G/T (rge) C/T? (90) C/T (85%) 
T = 29K .0113 .0097 . 0108 
3 20113 .0100 . 0112 
4 :0113 2 .0115 
5 .0113 .0103 .0116 


(b) Reese and Tucker [ 34] 

The results of thermal conductivity measurements on three sam- 
ples of polyethylene were used by Reese and Tucker to indirectly 
determine the specific heats of the samples. This led to a low pre- 
cision averaged set of esaet the temperature range 1 to 4.5°K. 
The comparison with the present work was made in the same manner at 


it was for Isaacs and Garland's work, and equation 4.1: was again used. 
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However in order to compare the values, a conversion was necessary. 
All of the results of Reese and Tucker were in terms of volume per- 
centages of crystallinities and the specific heats were in units of 
nj/OK*-en3. The values used below are not as reported, but have been 
converted using their values of the densities. 

The present investigation was not independent of the work of 
Reese and Tucker, and, in fact, was an outgrowth of it. Therefore 
the same types of polyethylene were used in both experiments.  How- 
ever slightly different densities (crystallinities) were found for 
the samples of the two investigations. This is explained by the dif- 
ferent molding processes, and cooling times, which is known to effect 
the degree of crystallinity [31]. 

The comparison is shown in Table V. The sample names are from 
the work of Reese and Tucker and should not be confuged with the same 
names used in the present work. The calculated crystallinities are 
shown in the table. 

These results are in agreement within the reported error given 
by Reese and Tucker (10%) except for Marlex I. However in view of the 
indirect methods of their measurements, it is not considered exces- 
sive. In addition they report a "dish shape" characteristic if the 
values of C/T? is plotted against temperature. This was proposed to 
be due to a lack of precision in the temperature and thermal conductance 
determinations at the two temperature extremes of their range. If it 
is assumed that this was the case at only the low end, then the points 
above 2°K will display one half of the "dish" and will appear as the 
curves in figure 3 between 2 and 4.5°K. This is believed to have 


been the case since at the lowest temperatures, the parameter from 
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which the heat capacities was deduced (the thermal equilibrium times) 
would be greatly increased due to the thermal resistance at the sam- 
ple-low temperature sink interface. In light of the above, the 


results are considered to be in good agreement. 


Table V 


Comparison of present results with those of Reese and 
Tucker (Values expressed in mj/9K -gm) 


T (°K) C/T? (Re T) CT ا‎ 


Marlex I (X = 79%) 


2 JOTA? -0118 
3 .0147 .0124 
4 . 0147 TUNES 
5 . 0147 -0129 


Marlex II (X = 62%) 


2 .0166 .0149 
3 . 0166 .0159 
+ .0166 -0167 
5 . 0166 .0168 


Low Density (X = 47%) 


2 . 0210 207775 
3 20210 s 0190 
^ .0210 .0200 
5 .0210 20207 


(c) Sochava [44] 
As mentioned in the introduction, the specific heat measure- 
ments reported by Sochava for the temperature range 17 to 60°K are on 


a sample of unknown density. Since it has been shown here that the 
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specific heat below 30°K is strongly density dependent, comparison 
is difficult. However, if Sochava's data are placed on a specific 
heat-crystallinity plot from the present work and are adjusted to 
fall on the straight lines, a crystallinity of 53% is determined for 
Sochava's sample. This is a reasonable value for commercial poly- 
ethylene, and it is felt that, as far as the comparison can be made, 
there is agreement between the two sets of data. 
(d) Dainton, Evans, Hoare, and Melia | | 

Dainton, et al, reported measurements on the specific heat 
of two samples of polyethylene for the temperature range 20 to 300°K. 
Their two samples were of very similar crystallinity (78 and 79%), 
and due to the manner in which the data were presented, comparison 
could be made only at 20 and 300K. Again the crystallinity was de- 
duced from the reported density (it was given for only one sample), 
but in this case it was identical with that given as determined by 
X-ray methods. To make the comparison predictions from the present 
results had to be adjusted to 78 or 79% crystallinities. This was 
accomplished by plotting the specific heat versus crystallinity for 
the three measured samples at 20 and 30°K, connecting the points with 
a straight line, and finding the specific heats corresponding to the 
reported crystallinities. Table VI shows the results of this predic- 
tion, along with the results of Dainton's et al. The subscripts 78% 


and 79% refer to the two samples. 
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Table VI 


Comparison of present results with those of Dainton, 
et al. (Values given in mj/9K-gm) 


Dainton, et al Present work 

C78% C79% C78% C7 91, 

T = 20°K 66.8 62.1 62.2 61.6 
30 147.0 146.2 146.7 146.4 


Considering that the 20°K results reported by Dainton, et al, are 
extrapolated from 22°K (the lowest temperature of actual measurements), 
the agreement is felt to be excellent. 

The Tarasov Model [ 47] 

The present results show clearly that a strong specific heat- 
crystallinity dependence is evident in polyethylene below 309K. 

Hence no single set of Tarasov parameters can be expected to predict 

the heat capacity in general, and therefore a choice was made to sep- 
arate the contributions of the two phases, amorphous and crystalline, 
and obtain two sets of thetas. But they are not completely indepen- 

dent, as will be discussed. 

The density dependence does not extend much beyond the tempera- 
ture range of this investigation. This was determined from the fact 
that the results of Sochava and Dainton, et al, are essentially the 
same above 509K, whereas they had differed by over 10% at 300K. Addi- 
tionally data from Dainton, et al, ona sample of crystallinity 58% 
above 90°K show that the specific heat for this sample is the same as 
for the more crystalline ones. (Above 110°K this is no longer true, 
but it is due to effects which are of no concern here). This com- 
pares well with the indications of figure 5 which is a plot of the 


relative difference between the amorphous and crystalline results 
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Figure 5, Relative difference between the 0% and 


100% extrapolated data versus temperature, 
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below 30°K. From the figure it is clear that the difference will be 
quite small for temperatures above about 60°K, if the curve is extra- 
polated above 30°K. 

The disappearance of the density dependence at the higher temper- 
atures has a rather natural explanation also on the basis of the 
Tarasov model. The three dimensional parameters associated with the 
two phases are small, which implies that the three dimensional modes 
will be nearly all filled at not very high temperatures. Thus the 
major contribution to the specific heat should be due to the one 
dimensional vibrations, which are assumed the same for both phases. 

The assumption that the one dimensional parameters are the same 
in both the amorphous and crystalline cases is based on the fact that 
these parameters are associated only with the forces along the chains 
which should be independent of the configuration the chain has assumed. 
It is also within the temperature range in which the density dependence 
disappears that the specific heat becomes linear in temperature, as is 
predicted by the Tarasov model, and it was from this temperature re- 
gion that the one dimensional thetas were determined. The linear 
portion of the specific heat-temperature curve yielded only an aver- 
age value of the thetas, and more information was needed to separate 
it into the longitudinal and transverse parts. 

The assumed value of 16559K for the longitudinal one dimensional 
theta was based on the measured maximum frequency associated with 
the C-C stretch of the polyethylene chain as determined by Raman 
spectrum [15]. Since this parameter had been taken to be related to 
the maximum acoustic frequency along the chain, the use of what appears 


to be an optical mode may not be valid, but this seemed physically 
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representative of what this mode should be. In any event it should be 
pointed out that the predicted specific heat on the extended Tarasov 
model is relatively insensitive to variations of the values of ہر‎ 


and Ôm as long as the relation 


سے 


a f 
4,2 L. ات‎ 21 


is satisfied. For example, if the ratio of Giz to Or is taken to 
be .6 (the same ratio as assumed for the three dimensional case), it 
changes the specific heat preduction at 20°K by only .5%. 

One further remark must be made on the one dimensional para- 
meters. On the basis of experimental data, obtained by neutron 
diffraction, the upper limit of the bond stretching acoustic mode 
observed in polyethylene was found to be 545 = 25 cm7 [34]. This 
corresponds to a characteristic temperature of 785 * 359K. If this 
value is assumed for the Ó,, it predicts specific heats which are 
much too high from about 129K on up in temperature. This is inter- 
peted not as anything wrong with the value of the maximum frequency 
noted above, but rather an indication that there is no real physical 
significance of the characteristic temperature from the continuum 
model. 

The three dimensional thetas for the amorphous and crystalline 
parts differ from one another in the manner in which they should be 
expected. They are related to the forces acting between the poly- 
ethylene chains, and, since the forces are most probably of the type 
which decrease rapidly with distance, the less dense and less closely 
packed chains of the amorphous state should have smaller thetas. The 


three dimensional parameter deduced from the experimental data support 
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this qualitative description. 

The calculated values of the specific heat based on the original 
Tarasov model are not in good agreement with the experimental re- 
sults. As can be seen in Appendix D, there is a 5-6% difference 
around 20°K, although in other temperature ranges the fit is better. 
The extended Tarasov model, which separated the longitudinal and 
transverse contributions, was used since it gave an approximately 
2-4% improvement over the entire range when compared with the basic 
model. However in the final analysis the Tarasov model should not 
be expected to explain the specific heat of polyethylene any more 
than the Debye theory can for other solids, since they contain inher- 
ent errofeous assumptions concerning the lack of dispersion. In view 
of this, the fact that the Tarasov model fits as well as it does 
makes it useful as a first approximation for the specific heats of 
fibrous solids such as polyethylene at low temperatures. 

The Stockmayer-Hecht Model [19,46 | 

In order to investigate the applicability of the model of Stock- 
mayer and Hecht as expanded by Genensky and Newell, the specific heat 
of the 100% crystallinity case was the only one considered, as the 
model was suggested for a fibrous crystalline solid. At the lowest 
temperatures, the model did predict the cubic temperature dependence, 
but the predicted temperature region for which the specific heat would 


be proportional to To 


was not found. However it is felt that this 
was an over simplification on the part of Genensky and Newell for, 
even if their analysis were correct, it would not be rigorously re- 


flected by this relation. The reason for this is that they predict 


in some temperature region the specific heat should be 
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where A and B are constants (the first term is the contribution of 
the interchain interactions and the second is the result of vibrations 
along the chains). This prediction does not give any appreciable 
temperature region in which the specific heat is proportional to a 
constant power of the temperature. 
To test the model for the temperatures at which the specific 
heat was investigated, a plot of c/T?/? versus the square root of the 
temperature should have a region where the slope is constant, if the 
model is correct. Figure 6 shows that this does not occur in the 
temperature range below 30°K, and, since the slope is negative at 
the high temperature end, it is doubted that it will occur. 
However, as is also shown in figure 6, plotting C/T versus T3? 


does show a linear region (22 to 279K) which implys 


TE Ce ae TS 


This could be interpeted as implying that their analysis is incorrect 
in the sense that they do not anticipate that the linear temperature 
contribution from the vibrations along the chain would occur at the 


5/2 relation holds for the interchain 


same temperatures that the T 
specific heat part. However, this cannot be reconciled with the model 
as presented, because the relative size of the two contributions is 
wrong. 

It is nót possible in practice to draw from experimental heat 
capacity data reliable inferences regarding the detailed form of the 


frequency distribution, as the heat capacity tends to smooth out the 


individual features of the frequency distribution. Hence the compar- 
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Figure 6, C as two functions of the temperature, 
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ison with the model of Stockmayer and Hecht cannot specifically deter- 
mine where their analysis is in error. On the basis of the apparent 
5⁄2 relation for the interchain specific heat contribution (if the 
interpetation of figure 6 is correct), it seems that the essence of 
the model is correct, but that a reexamination of the one dimensional 
part is in order. However other assumptions of theirs may be wrong. 
Since the polyethylene crystal appears to be lamellae, (although not 
directly observed in the bulk polyethylene) and the crystalline form 
is spherulitic, the structure may be such that the idealized tetra- 
gonal crystal, assumed by Stockmayer and Hecht, will not yield the 
correct predictions. Moreover since they had to make estimates of 
the relative sizes of the force constants, of which little 1s known, 
وہ‎ assumptions here could lead to different forms for the fre- 
quency spectrum, and thereby errors in the deduced heat capacity. In 
any event in its present form the specific heat predicted from this 
model is not completely applicable to polyethylene below 30°K. 
The Excess Heat Capacity 

In the absence of any acoustical measurements made at low temp- 
eratures for the heat capacity of polyethylene, the humping of the 
C/T? versus temperature curve, figure 4, was interpeted as an excess 
heat capacity. The appearance of this excess was not expected, but 
it is not unreasonable that it does appear when it is recognized that 
it is present in other amorphous solids [4,12, $]. The origin of the 
excess heat capacity ig unknown, but attributing it to the existence 
of non-acoustical low frequency modes seems reasonable. In particu- 
lar, the Rosenstock theory is appealing since it is known that free 


volume (holes) does exist in amorphous polyethylene [10,118]. However 
همه‎ 
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no reasonable conclusions can be reached on this or any other model 
which might explain the excess on the basis of the present results. 
The addition of the Einstein term to the amorphous heat capacity 
empirically fits the data quite well at the low temperatures, and no 
other correction is deemed necessary to account for the excess in an 
analytic expression for the specific heat. [In fact, it is felt that 
when low temperature acoustic measurements of the specific heat of 
polyethylene become available, the value for the extrapolated T = 8۴ 
point will agree with the following relation which is based on the 


present results; 


4.5 c E TA) TER 


From the size of the constant which multiplies the Einstein 
function, it is readily deduced that of all of the possible vibration 
modes of the solid, only .17% are taken up in the excess term. This 
can be compared with the number of these modes in vitreous germanium 
(.4%) and vitreous silica (2%), but no conclusions can be drawn [25,17]. 

Summary 

As has been shown, both the continuum, Tarasov, model and the 
lattice dynamical, Stockmayer-Hecht, model fail to completely predict 
the specific heat of polyethylene. This is not surprising as they 
are based on models which cannot be exact for differing reasons. The 
continuum model is at best a first approximation since it ignores the 
dispersion of the sound velocity which surely takes place in any real 
solid. Hence the spectrum of vibrational states is wrong. But, con- 
sidered as a simple two parameter approximation, the predictions are 


very good. The Stockmayer-Hecht model which attempts to account for 
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the dispersion fails for different reasons. Since this model is 
based on an assumption of the interparticle forces in various direc- 
tions, there are many parameters which determine the predictions. If 
some are chosen incorrectly, the whole model fails. On the basis of 
the present results no attempt was made to evaluate the parameters, 
since one of the basic predictions of the model did not materialize. 

The Einstein term adjustment for the excess heat capacity found 
in the amorphous polyethylene was adequate to explain the departure 
from the cubic temperature dependence of the specific heat at the 
lowest temperatures. But until such time that empirical evidence is 
brought forward to verify some model to explain it completely, it must 
remain as an a posteriori addition to explain the experimental re- 
sults. 

The hypothesis that the specific heat of polyethylene is linear 
in crystallinity was proven, it is felt, in the temperature range 
covered by this investigation, and extrapolation to the completely 
amorphous and completely crystalline cases 1s considered valid. It 
18 also believed that the crystallinity dependence has disappeared by 
50-60°K (considering only the temperatures below 110°K). 

The present low temperature specific heat measurements are more 
reliable than those given by the prior investigations (Isaacs and 
Garland, and Reese and Tucker), since the existence of the excess 
heat capacity is recognized. By extrapolating these results to O°K 
and by using the method of interpolating between crystallinities 
noted before, this work can be mated to that of Dainton, et al, to 
give a continuous specific heat-temperature coverage up to 1309K, 


where many other data become available. The resulting information 
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should prove valuable to workers investigating models of crystalline 
and non-crystalline fibrous solids into which nearly all polymers can 


be classed. 
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APPENDIX A- Calibration of the Thermometers 
Introduction 

Probably the most exacting aspect of this experiment was the 
determination of the absolute temperature scale. The difficulty 
stems from the fact that between 5 and 14°K there is no convenient 
temperature standard. Below 5°K the scale had been set by inter- 
national convention in terms of the vapor pressure of liquid helium, 
and between 14 and 25°K, the vapor pressure of liquid hydrogen is a 
good standard [8,48,51]. Additionally the internationally accepted 
platinum resistance thermometers can be used above 109K 177 Be- 
tween 5 and 14°K the temperature scale is normally established by 
means of a gas bulb thermometer. 

In this investigation the resistance thermometers were calibra- 
ted below 4.2°K by means of a vapor bulb thermometer, and from 4.2 
to 309K by means of a gas bulb thermometer. Liquid hydrogen vapor 
pressure thermometery was not attempted due to its inherent dangers: 
the facilities were not available to obviate them. Initially it was 
planned that the platinum resistance scale was to be compared with 
the gas bulb results, but was not since another comparison became 
available (the specific heat of copper, Appendix B), and it soon ap- 
peared that the higher accuracy in temperature range above 109K was 
not necessary for the purpose of this investigation. 

The theory, construction and experimental procedures concerned 
with the calibration of the resistance thermometers is the subject of 
the remainder of this appendix, together with an estimate of the 


error in the temperature determination. 
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Theory of Constant Volume Gas Bulb Thermometry 


The equation of state for an ideal gas is, 
Ve 
۸.1 Bae 


where ۸ and £ refer to the initial and final states of the system. 

If the volume is held constant and the initial conditions are known, 
then there exists a linear relation between the temperature of the 
gas and its pressure. A constant volume gas bulb thermometer could 
consist of a large volume, V, at the temperature to be measured, |, 
and a pressure sensing device at temperature, ly , with volume V con- 
nected by a tube of small volume yr’. Then equation A.l must be modi- 
fied to account for the fact that not all the gas is at the low 
temperature in the final condition; some is still in the gauge and 


connecting tube. To account for this equation A.l is replaced by, 


+= T, 
eq s s iis کا‎ asit (e) 
“eT 


where re is the pressure of the system when it is 1 and P is 
the pressure when the bulb is at temperature | . The integral ex- 
presses the fact that the temperature along the connecting tube 
varies from | to lí. However, it is convenient to assign a temper- 
ature Te such that, 
E Ux ole '(«) | 
t 


As 3 i ام‎ 
کرم‎ 


Solving equation A.2 for |, it is found that, 


A.U ar c. MEE ME 70 
: I rz t (r PEN) 
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/ 
Consider now the case where V~ Ô, that is assuming that the 
connecting tube's volume is quite small in comparison with that of 


the other components. Then equation A.4 becomes, 


سے 08 
300027 


Then fee, and if the bulb is at low temperature, most of the 


۸. بپ‎ x 


gas will be in bulb, and the pressure will be essentially linear in 
temperature. 

However, in an actual experiment the gas will not be ideal, as 
implied in the previous arguments. The equation of state in this 


case can be written, 


f CL |. 
ao )مء ع‎ + BO + = + ) 


4 
where the B; s are the virial coefficients, and L^ is the specific 
volume. If helium were used in the temperature range of the present 
é $ 
experiment وګ‎ is of 2 or 3 orders of magnitude smaller than B, » and 
© . / . 4 ° 
the term containing B, can be neglected [22]. Since Ë, is also small, 
KT. ; 
(can be set equal to "p in the term which corrects for the second 


virial coefficient, so that, 


a7 PVemR(T +t BP) 


f 


putting B =B, 


Using A.7 equation A.2 becomes, 


V u— un! ۷ v u- / 
A.8 —— — د‎ — m سس پر سوت‎ 





: I 
where the subscripts on the B$ refer to the temperature at which 
they are evaluated. 


Finally, consider the form that equation A.8 would take if the 
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initial conditions were varied somewhat, that is, by filiing the 
system with the gas when the temperature of the bulb is at a third 
temperature T, and therefore has volume W due to thermal contrac- 


tion. Under these circumstances, equation A.8 will be, 


Vw ص‎ a 1 
A.9 o 7 7,4 B, FP, Uk eT m 








LR. = s E 
P (ss + T, +G, P ٩ TBP) 


where the Po refers to the pressure in the system when the bulb is at 

temperature T, . (Notice that the Tc's on the two sides of the equa- 

tion are not the same since the lower temperature is different.) 
Construction Gas Bulb Thermometer 

As indicated in the comment following equation $€.5, temperature- 
pressure linearity can be obtained if the gauge volume to bulb volume 
ratio is made small. As a design criteria a was taken as 1/20, and 
the connecting tube's 6 U^ would be held small compared to the 
volume of the gauge. 

To obtain the desired temperature range, the bulb was initially 
cooled to 4.2°K by a liquid helium bath, and then heated electrically 
above this point. The bulb was made of copper, to insure uniformity 
of temperature throughout it, and the resistors were placed directly 
on the outer surface of the bulb. 

The pressure measuring device was a Wallace and Tierman gauge, 
type FA-145 (2 turn dial, 0-800 mm. Hg range). A .44 mm. inside 
diameter stainless steel capillary tube served as the connecting 
tube, and the gas used was helium. 


The gauge was calibrated against a mercury manometer (assumed 
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correct) both before and after calibration of the resistors. There 
appeared a linear deviation with pressure which was regarded as a 
systematic error in the gauge and a correction for it was made in the 
subsequent calculations. 

The volume of the gauge,™”, was determined by comparison with a 
known volume. The value was 


£0.10 +12 K0 PB) eus 


where ° is the pressure in mm. of Hg. (This correction term accounts 
for the fact that the gauge had a bellows type pressure sensing ele- 
ment). 

Since it was desired that V the volume of the bulb should be 
at least 230 cm3. A 1 1/2" copper pipe, 10" long, sealed at the ends 
produced a volume calculated to be 273.4 * .7 cm? at room temperature. 

Because of the size of the dewar, it was necessary to have the 
capillary 48" long. This resulted in a calculated volume for the 
connecting tube of only .19 cm3, well within the design criteria. 

As liquid helium temperature (approximately 4.2°K) would be the 
lowest temperature used in the calibration, some means must be de- 
vised to thermally isolate the bulb from the bath, once it was in 
equilibrium with the bath, in order to raise it above this point. 

For this reason the bulb was suspended in a vacuum container which 
was to be immersed in the bath. The bulb support was a 3/8" O.D. 
(.01" wall thickness) stainless steel tube. This tube also extended 
up and out of the cryostat and served as the pumping line. An elect- 
rical heater of about 1500 ohms, consisting of 12 feet of .002" man- 
ganin wire کہ‎ 78.3 ohms/ft.) was wound around the bulb and glued 


there with G. E. 7031 varnish. The heater was necessary to raise the 
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temperature and to compensate for the resultant heat leak through the 
support (and to a much lesser extent, the heat leak due to radiation 
and conduction by any residual gas). 

The resistors to be calibrated were inserted into small recep- 
tacles in the bulb. In order to avoid the relatively large electrical 
lead resistance, a three wire connection was made to each resistor. 
This, and a discussion of the low power A. C. bridge used, is given in 
the work by Reese and Tucker, and will not be presented here [ 34 ]. 

A schematic diagram of the completed system in shown in figure 7. 

Gas Bulb Calibration Procedures 

The gas bulb was filled with gaseous helium to about 790 mm. of 
Hg when the system was at liquid air temperature. The vacuum spaces 
were then pumpéd to approximately 1077 Torr and the liquid helium 
space was filled. 

The data accrual was quite simple. The procedure was to heat 
for a period of about 30 seconds, then turn off the heating power, and 
then observe the resistances come to equilibrium, which indicated 
temperature equilibrium. The pressures and resistance values were 
then recorded. The tendency for the temperature to drift downward, 
due to the heat leak, was compensated for by the heater. 

The observations were taken over a pressure range of about 40- 
300 mm. of Hg, corresponding to temperatures of from 4 to 30 degrees. 
They were taken twice over the entire range, at about 2-4 mm. steps, 
to check reproductibility, and a third time over the lowest portion 
of the range where the resistors were the most sensitive.  Seventy- 
three points were taken altogether. 


After completion of this procedure, the bulb was allowed to come 
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a, Vapor bulb 


b, Resistors to be calibrated 


Figure 7, Schematic of gas bulb thermometer, 


to complete equilibrium with the liquid helium bath. Using the vapor 
pressure, the bath temperature was determined. This would later 
serve as a check on the accuracy of the calibrating procedures. The 
helium was then evaporated from the dewar and the space filled with 
liquid nitrogen (99.9% pure). The pressure in the gas bulb together 
with the known temperature of the nitrogen bath gave values of f 
and T, which could then be used in equation A.9, and the pressures 
could be converted into temperatures. 

Equation A.9 is here rewritten with a slight amplification of 


the terms. 


fu Y» ۳ U Fu! 7 
0ھ‎ Tu t BuPy ° Tv+4+By Ted eu ` 
PV Pu 2 





aap ۔‎ 7+8۶۰ CUA 

where f is the pressure in the system when the bulb is at liquid 
nitrogen temperature, 7, is the liquid nitrogen temperature, V is 

the volume of the bulb at temperatures less than 309K, and Vw is the 
volume of the bulb at Jy. The subscript ¥ refers to ambient condi- 
tions in the laboratory, and the other symbols have the same meaning 
as before. Using published values for the thermal contraction of cop- 
per at the various temperatures, it is found that V=9270.81.7 cm? and 
E 273.01 .7 cm? [26]. For an estimate of the Due » a value of 145°K 
was assumed for ns and 1089K for 7.2. These were obtained by assuming 
a linear temperature distribution along the capillary. This is suf- 


ficiently accurate as this part of the system is quite small. Values 


of the second virial coefficients were taken from Keeson;[ 22] 
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B = (.2453 - 6.2009/T) °K/mm of Hg‏ 1ھ 
B, = .000163‏ 
Gy = .000162‏ 
e 2 .000179‏ 
c 7.000186‏ 


Using the above in equation A.10 leads to an equation quadratic in 


Solving that for T yields, 


۸.2 T= --« 4l 4 - 4a 
where a = us 
- (€: - «)P 
ce late + cag + oz) - FE" = 
olen جب کے‎ Bz G.2009 , 


Equation A.12 was programmed for the CDC 1604 digital computer. 

The procedure discussed above gave only a pressure-resistance 
relationship until equation A.12 was solved. To obtain the desired 
temperature-resistance dependence, all the data was best fitted, in 
the least squares sense, by the computer to the following type of 
equation for each resistor; 

7 m 
A.13 T= £ A: (Ly. 1 ° 
Axl 
where the Ais and the best fitted ^ (integers between 1 and 7) 
were determined. The values are presented at the end of this appen- 


dix, and the RMS deviation from this fitting is calculated to be .23%. 
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Calibration of the Thermometers below 4.2°K 

For temperature determination from 1 to 5.2°K, the International 
Commission on Weights and Measures ad@&pted the "1958 He Vapor Pres- 
sure Scale of Temperatures" as the international standard [ 8,48]. 
Tables based on this scale give a conversion from vapor pressures to 
temperatures, hence the temperature can be ascertained from a reading 
of the vapor pressure alone. 

For this lower temperature determination, it was decided to use 
a vapor bulb thermometer, rather than merely measuring the pressure 
above the bath since the latter procedure involves making corrections 
and assumptions of questionable validity. 

Vapor bulb thermometery is performed by placing a container into 
the helium bath at the liquid level at which the temperature is de- 
sired to be known, and filling the container with gaseous helium 
until the gas condenses which indicates equilibrium between the 
liquid and gas phases at this pressure 2سا‎ The pressure is then 
observed, and tables give the temperature directly. 

For the actual construction of the bulb a small copper cylinder 


3 was connected by a 1/8' stainless 


with a volume of approximately 2 cm 
tube to a mercury manometer and filling manifold. The bulb was 
mounted directly on top of the vacuum container, and connecting tube 
was brought outside of the cryostat enclosed in an evacuated tube (to 
prevent gas oscillations in the system [42]. See figure 7. 

For the calibration, exchange gas was introduced into the vacuum 
space where the resistors were located, and the pressure above the 


bath was lowered in steps. The manometer readings were taken and 


the corresponding resistances were recorded. In all 32 points were 
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taken in the temperature range of from 2.4 to 4.2°K. 

After conversion of the pressures into temperatures, the resis- 
tance values were related to the temperature by equation A.13. 

Recalibration of the Thermometers 

Inasmuch as the comparison experiment (Appendix B) indicated that 
a systematic error might be present in the temperature scale, the 
scale was reestablished. Also, in the initial calibrations, the vapor 
pressure and the gas bulb temperature scales (below and above 4.2°K) 
were independently determined with the result that there existed a 
slight discrepancy at their common temperatures. (Approximately .2% 
difference). Since it was desired to have a continuous scale over the 
entire range, the calibration was repeated insuring that the tempera- 
ture overlap was consistant with both types of temperature determina- 
tion. This was accomplished by using the liquid helium temperatures, 
as determined by vapor pressures, as the fixed points for the gas bulb 
calibration rather than the liquid nitrogen point. As a check as to 
the precision of this method, the predicted liquid nitrogen tempera- 
ture was only .39K different from the assumed value of 77.359K, a 
difference of only .23h. 

The equipment and procedures were the same as outlined earlier. 


A total of 99 points were taken over the entire range, varying the 


i 
رڈ‎ 
F. 


intervals so that a uniform coverage was obtained. The gauge volume 
was again experimentally determined, and was within 3% of the original 
measurement. The gauge was likewise recalibrated against a mercury 
manometer twice. 

Ihe resistance values were related to the temperature by use of 


equation A.13. 
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Accuracy of the Calibration 

In the temperature range below 4.2°K where the vapor bulb ther- 
mometer was employed, the only significant source of error in the 
temperature determination was the reading of the manometer. It could 
be read to within .001 cm of Hg with reproducibility of .003 cm. of 
Hg. This results in a maximum uncertainty in the temperature for this 
range of 3 millidegrees. 

For the gas bulb calibration, which covered the temperature range 
of from 4.2 to 30°K, the following effects were considered to ascer- 
tain an estimate of the error; 

(a) The fixed temperature point. 

The calibration over the entire range is very sensitive to 
any errors in the determination of the fixed temperature point which 
essentially sets a value for the right hand side of equation A.12. To 
determine this parameter, six independent observations of the tempera- 
ture were made, by vapor bulb methods, at the liquid helium bath tem- 
perature, and the pressure associated with each recorded. The fixed 
point was taken to be the average of these data. The resulting error 
in the temperature determination over the entire range due to this 
was estimated to be at most .003 degrees. 

(b) Errors in the measurement of volumes. 

Since the bulb's volume was much greater than the other vol- 
umes in the system, the temperature is essentially linear in the 
pressure. This tends to make any error in the measurements of the 
volumes self compensating, and their net contribution to the error in 


the calibration is negligible. 
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(c) Thermalmolecular correction. 

Following Roberts and Sydoriak, the thermalmolecular cor- 
rection is .004° at the low end of the temperature seale (4.2°K), 
and negligible above 7°K [71 This systematic error could have been 
corrected for, but was not in view of the existing larger errors which 
will be discussed below. 

(d) Pressure measurement. 

The pressure measuring gauge readings were compared with those 
of a mercury manometer both before and after the resistor calibration. 
This comparison was made stepwise going up in pressures, since this 
was the way data were taken, to include any hysteresis of the gauge's 
bellows as part of the correction. The gauge was determined to have 
a linear error as a function of pressure, that is, the deviation was 
consistently high at the higher pressures, and low at the lower ones. 
Ihe data were corrected for the presence of this systematic error. 

The statistical uncertainty of this correction is entirely compatible 

with the precision to which the gauge could be read, about .1 mm. Due 
to the uncertainty in the pressure determination, a possible error of 

.019 exists throughout the temperature range. 

Other sources of error, such as those considered by Beattie, and 
those due to the uncertainty of the values of the virial coefficients, 
accuracy of the mercury manometer, etc., were considered to be insig- 
nificant compared with the ones discussed above [5]. 

Once the temperature scale was established, there exists the 
necessity of relating the resistance values of the resistors with the 
temperatures. This introduces another possible reading error which 


amounts to .1% over the entire range. However it is felt that due to 


the insensitivity of the resistors at the higher temperatures, a more 
realistic value of .3% should be assigned there. After the resistance- 
temperature relation had been established by equation A.13, it was pos- 
sible to assign the maximum temperature differences these errors would 
introduce, .00159 at 2.5 and .039 at the highest temperatures. It is 
also noted that to use the resistors as thermometer, another reading 
error of this size could exist. 

Sunmarizing all of the above, and recognizing that the fitted 
relation between the temperature and the sample thermometer resis- 
tance values contains an uncertainty, it is felt that the temperature 
determination by use of the resistors was accurate to within .5% of 
the absolute temperature. 

As an estimate of the reproducibility and precision of the cali- 
bration, the RMS deviation of the actual points from the fitted 
resistance-temperature was used. This amounted to .23%, and appeared 
to be randomized. Hence the statistical scatter of the calibration 
points from the fitted relation was well within the estimated errors 


of the procedure. 
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Parameters for the temperature-resistance fitted relation 


T> É Ai (hoyu R). 
Thermometer 
Shield Shield Sample 
(4. 2-30°K) (2.4-4, 20K) (2.4-30°R) 
n 6 6 i 

Ay .81986x107! ~.11574x1071 -.20787x109 
A2 . 6500X1072 .66904x1071 . 32888x10! 
A3 -.50749x10! .25718x10-l -.10029x10¢ 
Ay .17961x10? .10939x10? -.3639ux10? 
As .14366x10? -.27080x10! .25491x103 
6ھ‎ ~.15249x103 -.69149x10! -.45361x10? 
A7 .19880x10? 62ء‎ . 36772X103 
% RMS 

Deviation . 20 . 00 23 
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APPENDIX B = Comparison Experiment 

There has been considerable discussion about the need for, and 
choice of, a suitable calorimetric standard in the temperature range 
below 20°K [53]. Because of the large number of specific heat deter- 
minations for copper, its available purity, and other reasons, there 
is much support for its choice as such a standard. At the same time 
that this experiment was being performed, Martin reported the specific 
heat of a cast 99.999755 pure copper sample in the temperature range of 
this experiment [30]. He asserts that his results are accurate to 
better than 1%, which implys a temperature determination to within 
. 35% accuracy. 

To ascertain an estimate on the accuracy of the temperature 
scale used in the present experiment, it was decided to duplicate 
Martin's measurements. A Single crystal copper sample of 99.9995% 
purity was obtained (Semi-Elements, Inc.), and was etched with nitric 
acid to a size which was compatible with the calorimeter. The re- 
sulting sample was a 24" by 9/16" rod with a mass of 83.181 grams. 
Heat capacity measurements were performed on this sample, in the same 
manner as for the polyethylene samples, over the temperature range 
US to 309K. 

In performing these measurements, the only experimental procedu- 
ral differences from the polyethylene runs were that the amount of 
heat introduced for the temperature steps was considerably less, and 
the equilibrium times were much shorter, due to the differing heat 
capacity and thermal conductance. Three complete sets of measure- 
ments were obtained over the aforementioned temperature range, re- 


sulting in a total of 119 data points. 
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After smoothing the data by least squares fitting it to a seven 
parameter equation (equation 3.2) and correcting for the heat capacity 
of the addenda, the results were compared with those of Martin. They 
matched fairly well at the low temperature end, this work being about 
2% lower, and excellent agreement was obtained at the highest tempera- 
tures where no difference was found. However Martin's results were 
over 5% higher at the midpoint, 17°K. The differences along the scale 
were not random but showed a clear divergence up to a maximum, and then 
converged to zero at the highest temperatures. This type of deviation 
can most easily be explained as some sort of systematic error in the 
determination of either temperature scale, or possibly a drift of the 
resistance values of the thermometers. Martin discounts this latter 
possibility in his work by repeated checks against both vapor bulb and 
gas bulb thermometery. In the present work nearly every time the sys- 
tem was at equilibrium with the helium bath, the thermometers were 
checked against the vapor pressure temperature of the bath. On the 
basis of this, it is felt that there was no resistance drift of the 
thermometers used here at least at the low temperatures. 

The size of the maximum difference in the two sets of results 
is quite large, but it should be pointed out that at 179K an error 
in the temperature determination of only 1.3% (.29) would explain 
the difference. This value of error is outside of the estimated 
error in the temperature scale used here, however, and so it was de- 
cided that a recalibration of the thermometers would be performed at 
the conclusion of the specific heat measurements on polyethylene 
(the copper data were taken about half way through the polyethylene 


work). It was felt that with the low temperature techniques acquired 
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by this experimenter that a more precise calibration might be achieved. 
The recalibration is discussed in Appendix A. 

After the recalibration was performed, the heat capacity data 
for copper were recalculated, and again compared with the results of 
Martin. The same relative deviation persisted, however it was approx- 
imately halved at the 179K point (2.6% low). At the highest temper- 
atures investigated, 309K, the results were again the same as Martin's, 
but the same 29; difference was still present at 59K. 

The new calibration was decided to be the more accurate and was 
used based on the following considerations. First the agreement with 
Martin's copper data was better. Next the recalibration was per- 
formed after this experimenter had considerably more experience with 
the low temperature techniques, and was able to look for small sys- 
tematic errors which may have been overlooked the first time. Lastly 
this calibration was continuous over the range 2.4 to 309K which was 
highly desirable (see Appendix A). 

At the conclusion of this work, Alhers published specific heat 
data on copper over essentially the same temperature range as that 
under consideration [ 1]. His work compares very well with that of 
Martin, but has the same relativei-^ deviation from this work. For 
comparison, the present copper specific heat data are presented be- 
low along with that of both Martin and Alhers. Additionally the 
parameters for the fitted curve equation (3.2) are given, along with 


its o/o RMS deviation from the actual data points, (.69%). 
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Parameters for the specific heat-temperature fitted relation 


7 t 
A- | 
C = pa A. T (Gin m j/9K-gm). For T between 2.4 and 30°K. 


As d 


Ay ~.12463X10° 
A» .96209x107! 
A3 -.23948x107! 
Ay -41835X10~2 
As -. 26648X1023. 
Ag .10385x10-^ 
Ay -.14265X1076 
% RMS .69 
Deviation 


Molecular weight of copper 63.57 gms/mole 
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APPENDIX C - Raw Data 
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حا اک نر تسس بب بو U 7 CU‏ ا ا یہر — ہے E‏ 
10.61803 9.990 , 

10.265 11.22138 

ج 3 0 ےر ہے ےج ےر ہہ 
12.30049 10.604 

10.834 12.99738 

aia روي قو ص0‎ ~~" Sao 
11.148 15.03310 

11.339 0+7 

pM — tO JA ۴نل‎ 
11.644 15.82448 

16.29977 11.786 
7 اوک وب سج سٹو i  ———‏ 
17.62685 12.162 

12.198 17.79657 

مک — روو و Kum 0 o9‏ 
20.48844 12.859 

13.239 221 

aa”‏ دک ہر ola‏ ہر یو ہہ 
23.74350 135 

13.856 24.79318 

En e a 77 ج په‎ 
14.343 . 95 

15.563 21٦ 

O 55)‏ —— ص ہے رہہ 
31.75482 14.973 

15.371 31.86838 

pu =. ee 
16.320 37.06861 

17.374 42.64875 
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MARLEX 1 


oe 


SPECIFIC HEATS IN MJ/OK/GM 


TEMPERA "XE o0 


48.82613 ۱ 18.430 
ITU. HO OTO‏ 
(OT‏ 15 18.875 
54.50185 19.274 
سک کک کے ا سس سم سوہ وکا مک ~ چو کے ٦.9‏ 
57.71516 19.799 
59.66975 20.105 
OD PTS‏ 
65.13815 20.884 
68.59212 5 و21 
ae — — 2l.55 [511‏ 
71.80872 21.919 
77.27470 22.568 
O= FF‏ مت سا سو کہ کک کا کک کک ار ا 
81.65808 22.828 
85.01005 23.499 
91.55677 24.186 
9585445 24.980 
Eg: DI a ONTO ————‏ 
100.1683530 2292 
106.23490 26.191 
—— سس a‏ )يئ )اس —290.5T8———‏ ——— <_< 
1167912352 27.159 
121.90134 252 
OOOO‏ وا و ہن مک کک nis a‏ 
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MARLEX 2 


SPECIFIC HEATS IN MJ/OK/GM 


——————————————  . س‎ A 


ao soem — — — TEMPERATURE یب سک‎ 


1.10836 9و یا 3ے ہا 
es SOS‏ .و رہ ہے ہہ ہک ee‏ 
WAT °` 1.15943‏ 
1.26999 4.538 
«eee cc‏ 540 
1.31760 4.568 
1.41884 4.703 
٢ 7‏ ) حتتت لل ې )0 Een e‏ 
1.50724 4.762 
1.60353 4.881 
a Z. — — 7 S — — — mel‏ 
1.70558 4.967 
1.82800 5.092 
ccu-- —5c OM — — — iss r‏ — 
f 5.187 1.96224‏ 
2.05350 5.301 
SO) Hs 0‏ 
2.22102 5.484 
2.25201 5.112 
a TT Og — 2 UNS) |‏ کے کک .مہ 
2.51030 5.670 
2.53578 5.645 
Gg SP O‏ فی مک یک ا اکا رر ا ےچ ہہ 
2.78156 5.859 
2.89658 5.897 
eG O O0 DS‏ 
۳ ولد 6.197 
3.31436 6.165 
5 سی کک کچ و6 
b454.. 3.75170‏ 
3.93908 6.595 
57 62118 
4.29967 6.770 
4.57285 6.952 
نت UI) Seo‏ 
4.93654 7.118 
5.26753 7.315 
ROCIO‏ 7 330 7 
5.69568 7.589 
5.90489 7.642 
eee o SS‏ — ر ي و 
6.505155 875 
6.69330 7.990 
oe — —— tia‏ = 
7.47455 8.338 
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MARLEX 2 





ہے  —‏ — ےہ دپ ےج وم ا سے 


SPECIFIC HEATS IN MJ/OK/GM 


TEMPERATURE  — -—— ct — 


6 0000اک حسسمججججججچ ےت وکا 150 


18.345 


97 


55295 


8.612 8. 36381 
M  GeofS 89 ^ 
8.720 E2132 
9.2164 9.63058 
————_ —9-5351 (162 03 U 
9-447 10.62664 
9.7111 11.613580 
9.907 oS —— —— —— — 
10.000 12.28690 
10.314 13.08954 
105 So su 
16.534 14.10426 
10.860 15.035178 
mot — ——F5*515062 
11.098 16.19879 
11.591 16.97032 
Vtr69g- — — 859 
11.694 18.33617 
171.938 19.13534 
126195 PO ONUS — — — 
12.315 20.92201 
12.495 21.41557 
12.797 22 بال ا72‎ 
12.9651 23 0 0 
13.066 25.89595 
ہے ےت‎ — 
13.701 26.843515 
13.660 26.951060 
1.۰1 با جا‎ 2382655 — ——— 
14.353 50.01182 
14.522 39.75287 
Mo: o SOU سب‎ 
15.0 30 33.22668 
15.236 34.73176 
15.738- 37.158855 
15.784 37.24379 
16.060 39.29161 
16.573 41248706 
16.596 u2.02919 
16.911 4.06952 
17.397 46.280928 
17.448 46.91189 
JT 7ھ‎ 49.84552 





MARLEX 2 





SPECIFIC HEATS IN MJ/OK/GM 


E يو‎ — — TMNPERATUTE — — CC 


— "amsa ME یمیس‎ —€—— 


18.832 55.59494 
coco C B.789 OL 
19.311 58.92028 
19.688 62.0974] 
وہ جک و رجہ‎ Ghost ON وو لاج‎ O> 
20.324 66222665 
202554 68.59343 
Ecco a یچچ‎ ٠ : ےو ویوووی‪ وڈ ڈو وہ‎ 77 
21.507 73.83893 
21.511 75.11742 
= e MELSUSO  - ٢م‎ ge 
22.729 0289951 
22.047 84e 30054 
co] 8) 20 ar 
23.871 92.26315 
23.844 96.20608 
ieee -  @¢4e200° ٢-0 کے‎ 
24.965 192.41275 
25.217 104.56977 
Z35469 10232250 
26.118 112.90491 
26.716 117.49972 
2051917 ~ T ba 
27.999 1512556095 
28.399 131.63365 
28,601 29م د1[ ۔.۔.‎ 
29.244 137.68534 
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LOW DENSITY 


س ا 


SPECIFIC HEATS IN MJ/OK/GM 


TEMPERKTUTE — م‎ C OS 





a 
—————————— 


— — ——oAy v لا‎ 


a M a هم‎ 
. 


]255 2.521 
MEME — 2.553 © «29006‏ 
31867 . 2.59 
335245 0 ]2.62 
amm CC OO e345‏ 
37404 © 2.695 
e 38825‏ 2.731 
o o oOo Sg rl, a TAT‏ 
042907 2.816 
Oe 7 7 7 uq‏ 2 — کک کا 
50078« 2.941 
51675« 2.978 
"ANUS I» NU)‏ 
moo‏ $5 3.097 
059378 3.088 
6ء یکو ہا ...سس سس شش — << 
)0772 5 
269384 3.266 
BT‏ — )لو اين 
81595 3449 
82114 ))باباء د 
وو ا 
e 96623‏ 5.648 
۱١٩ا)ه.‏ 3.657 
— جحد" ١‏ د مک a‏ 5 3 
1.14367 3.851 
1.14469 3.863 
1ر چا ہ چو Og e‏ 
1.33820 | 4.061 
1.35516 4.064 
t.2T8——————17457T932 —‏ —— — — ص 
1.61775 4.328 
1.62255 4.318 
s‏ ونو د يڅه 
1.81921 ہا 7 بجاے با 
1.81792 4.498 
D d coe — cI DES 58 —-‏ 
1.86155 4.508 
24035621 4.642 
es — —1.682——————————— —2.00809‏ 
742 2 4.706 
2.30063 4.831 
سور سو = کک کک ا ہی مه 
26.46734 4.932 


— — 


29 


Loe 


SPECIFIC HEATS IN MJ/OK/GM 
"=" TEMPERATURE چپ‎ 5 ٣٠٤٥٠ 





F - 


— ——T.———— —— naa — د‎ 


100 





روسمس ——————€ 


5.034 | 2.63293 
LE —— 0 o S$5.lió © e ١ 
5.1506 2.79630 
5.189 2.88276 
= 5.205 —— oou c c 
5.360 3.19051 
5.382 3.21864 
LEN DE o DUDU — "vuol 
5.470 3.36955 
5.615 3.64384 ` 
Pee ee bau s mn m 
5.712 ` 31.831 
5.713 . 3.85035 
| ca i قحلا ححست وويم ون‎ | 
. EM 5.916 4.27282 
5.980 4. 32430 
MEE. mes 72 4.45900 
6.209 4.91476 
6.235 4.95730 
5:27 ځ‎ EEL E ST 
6.303 5.10324 
6.550 5.69680 
6:5715 S755 
6.592 5.78959 
6.644 6.00095 
bau” OTTO 
6.928 6.63617 
6.922 6.65371 
lu EE ————— — 17.268 "^ —"'Tra52uT2 
| 7.263 7.5614) 
7.287 7.61844 
T6235 -~ b] 
x 7.669 8.71615 
| ۱ 8.024 9.80153 
Wa o ون > سسس ص‎ LC Oe 
8.505 11.061 
8 با ء‎ 6 11.21126 
رہ ڈ ہا ری وہ ہے رہ‎ IT ey S a o 
8.815 12186119 
| 8.899 12.75976 
i Se es cages Son CC aar 
9,361 14.4457] 
0 See. ee ae 
9.661 15.61894 


سو — مر بس — — 


LOW DENSITY 


SPECIFIC HEATS PRE MJ/OK/6M 


iii” TE TEMPERATUXE 8 «©  G- 








ج 


101 


| 9.884 16.52152 

a —s».89g o T65?95 
10.127 : 1759983 
10.464 18.91642 
TO. Se 18.95116 D pgs 
10.639 19.73098 
11.038 214275 
OG 2109968" 
11.195 22.29614 
11.631 24.39704 
11:669 24.42918 
11.788 25.18032 
12.229 27. 30692 
IN me 27.65680 
12 10 28.33392 
12.845 50.68352 
12.901 30.86788 
13.063 31.5591 
ES 34.28889 
13.576 Soo. — 
13.767 35.76025 
14.245 38.38456 
U Teis 38.66656— 
14.528 39.96585 
5 22 42.96486 
15.098 با‎ 373253 
15.340 44. 96633 
15.862 48.09940 
15.940 L8. 475]] > 
16.220 50.30252 
16.758 53.8299u 
16.845 5,761) سے هم-5‎ 
NE 56.28382 
17.700 59.96677 
WIT 65 60270527 
18.713 67.33329 
oum 67.69524 
19.002 69.52506 
19.788 75.32322 
19.860 ھ3۰٦۲‎ 
19.986 76.79121 
20.901 84.22005 
20.958 84.50975 
21.055 85 ہج وہ‎ 
2 DU f 92.92935 


LOW DENSITY. 
SPECIFIC HEATS IN MJ/OK/GM 
lac TEPE eG 


22.118 94.185537 
22.118 94.42999 
25:228 10525555! 
23.3556 103.53852 
Zoveil —— ^ "AOUmnd5DUS O © 
24.385 115295511 
24.352 | 115.23389 
24.560 116.08967 
25.490 124.87124 
۲ء256‎ 125.00688 
25/6 125.969»T 
26.735 133.83623 
26.654 134.543406 
Exc  —2060.992^ - - vu» 
216869 140.81310 
27.939 143.93052 
یتہج تا‎ «26.215 o a ESOS 
| 29.122 154.61557 
| ۱ 29.283 156.00964 
29.038 61.87336 


—KL—sv"—1aƏ‏ ےسیو سے 
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وسين نم Y —— — - — . —— a‏ ےب نے 


CRYSTALINE SAMPLE 


THE 100 0/0 


CAPACITIES IN MJ/UK/GM 


HEAT 


TARASOV C 


C ACTUAL 


TEMPERATURE 


DAP BOOM ODT Ne POR DF ONO ۰۲۷( یہی یی ں۔‎ IT JM یىی تہ‎ O OA C O N< 3 =— 2ه ہے ہی ایہم ری ری‎ ۲0 FN ON ath 

o 9 o e 6 @‏ د چج 9 9 ده هده 9 9 0 99 9 0 9 .9 9 0 9 9 9 6 9 *9 P 09 9 9 9 o@ 6 6 6 6 o‏ ° 6 هه 909 9 هه 90 @ ده 0 9 © © 

C06 ١ | CA bee NNN ONS =+ + + + + + + + > 31 y + + yrOOrOOCNir N c | MUNIN O 
|| | |١ || 5959014507 || l l ١ | | | | ١ | ١ | | ١ | ١ ٠ ٠ ا١1۹ ہم‎ 


MOM TMH OCC FrOCIMNOR MOMMY DW = Mr hOM MOAIMMIN =H MAIN Ft ONE Fre MONI RH OOM 
Oe 0) ما -— ۸ يم مه همر سم,ن‎ ۷( IC COO سو الہ ن سوج ہہ ری ن۔ ن)۔ رح یہ‎ TIM ین‎ MOAT FOTONO = LAMM OOO ha re 
(یا ۳۷ یح سم‎ ٣- ن٠٥ رح ۶م م سم ى۷۲‎ OO [a MA O — ICG IO هټ‎ MAIN DN AACR OD wae ONE CIO ARIANA O truy r 
e 9 همو ه ه ه مج ه ه 9 90 09 9 9 09 یم 009 0 09 ه هه د 09 .9 9 090 هه 09 0 0 09 9 0690 0690 90 0 ه٥ 090 09 9 90 090 ه ه‎ 
~ere“ (NIM TF run CR 0 Cor rou - C سم‎ م۸٣‎ 07 O rr: U سس همپ سم(0((‎ —+ I— — + سم,نین‎ ١ دن رم‎ A DONO + Or 

ېم + O: C G Ç O°‏ م000 OOO EL‏ +۔ (۳٠٩٣۷ ( HF >ü‏ ټم ڼم) زح لپ سر سم سم سم سم سم 

r= r= N OO ram 


pr 


- 


THAIN HN TNC SNE cic uxo orum o) Xo NIC CCOJA Ur pE-uit-ooOormgqcoel o9. —-o0zroocogo 
Oe PO سم‎ Or Of c7 PO cf TMH ODM OO FH TOM OUD سو حم پہٰ‎ ce rer 0 م 676 )هم من سم‎ IT بت ه‌ت ه ی٣ ډ يم )۱م‎ 
(ا(۴ لہ سم‎ ٣- ۳ عم‎ re NID COPE OOOO OUI NII ONE) SOR م‎ OON Oraa [Mm Ora TONG DH GO 
€ 99 0€ 90 90 € 9 06e 0 6 6 6 6 6 90 0 9$ 0 0 € 9 09 6 هه‎ O O 90 09 0 6 6 6 9 00 9 9 0 00.0 0 6 6 6 6 9 6 دو ه‎ > 
س, سم‎ NNO +۶ ی‎ Of- CC CD 0 رہ ں ۶ ہ۷۸‎ D 0 ++ OQ e- zEP-OS CU OOCN UY OO CNUOY O 0-0 05 zr OS OP- CONI GOGON CN 
ecce oo qJeogm N zr xr -unwunuo-o-Or-h-0200000 OC2 Orngoo —+ Os 

: rare N) O Ge 

p= 


NIN GIP QANDNOMN OWNS NOUN ANDO NOMNADNOQMNAN SinDnNn OMOMOWONTVMODNTMNONODCOO 
690900606069 090620999 نه‎ O 0 هو‎ o O وه د هن هد هن نه ° دن نه م ن هن @ ده ده @ 0 و دنه ےم 0 وو 0 @ دن ه نه 0 هن هه د 0 نه هه‎ 
NIN SF TUDO OR یی زی رک مخ‎ l ت‎ —cjJ(NrOrO iy م ټ ۷۵ ی+-‎ DM OO GO emen NGGON نن هه ټم ټی يیایا:‎ 

e onm Kan — r— r r— r— r— r— سم سر سر مس مس مر سو‎ o (N NJCN TICN CIN اہ لہ یہ یہ ی0[)) یہ یہ‎ OUO 


r—r— (NJ 
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THE MARLEX ONE SAMPLE 


APACITIES IN MJ/OK/GM 


po 
تیه‎ 


HEAT 


TARASOV C 


C ACTUAL 


TEMPERATURE 


NON Mm + — COCCO DOrn— nO OC I IO cJu- De یی‎ ۱00۸6 I N i — O P— —-+r - O د له خا پړم .کپ ڼ ده‎ OCI ONMO 

هوه یه هو ه هد و هو ه 09 و هو ه هو هد هو د و 09 ه o 0 ٥‏ 0 ه 0 و ه هد و هو ه 09 o o‏ 09 0 هو مه 09 o‏ دو 09 ه 09 0 ده ه م نه هن ٥‏ ه هن 

I ١ ١ | | | | eK JNNMMM + + 9 (١ ١(۱/()١(ا١( کا‎ (۱ ١( 9ے‎ + + ٣ ٥۳(۱۷(۷۴( ابم ورسم سم 141 سمسمیہ لح‎ 
ll l lil 1 l jl l l j l f lal l l ll l j jp l NGNE نے‎ 


ON OM OLN AF N IOOC ONN.‏ سر ہے ے زی ہی یی CC OO‏ ین یملںم ےم NINA AN eR Sm OWN OM‏ ۶رز ںن زی می نب ۲۴ سم 
هف ېي خنں ۳۱ہ یہ(٢۷۱۳ایہ)‏ ئک تن )یم نلم (۱۱۲ہل)) R= DOM + ON AE NI POND f CO ze e MO OR SDE FIND OR‏ 
COR OR AN OH ON ODA AOR OODOR-ONDENOOAIM OM r= OF ON AMNOMO‏ سم +- 0۰ +۲ — هم += r‏ — 
هم O‏ 6 50 6 6 6 6 6 6 6 9 6 6 9 9 6 6 6 6 6 6 6 6 9 9 ه٥‏ هه 9 9 0 9 9 9 090 0 9 9 9 6 6 6 6 6 6 6 6 6 6 6 6 6 e o‏ 
ON GO AM — 27 OOOO FOOr-O‏ ال" ری ص) MH‏ رآ ال۱ حن NEI Ot LD OR OT CNM Orme N OOO =I O‏ سم 

rara amma ONO AISI IFO + + ۶ می‎ (0 -O O= == O O O O: O ی سم تت‎ FN 

O" —‏ لپ لح سس p‏ سس سے 


gaman 


~ 


Chi rn OUINIC COO HOMO ONANNOMA-OGOO‏ رت COLO y -- At OOR MME DR IO e e NH)‏ ۲۳: ۱۷ہ ۸ہ یم 
NO OOO‏ ۳ب ۸م سیہ هی سم 0۰۲٢۷‏ ۶ے OM NO TFN‏ 0م ی م CJO P un zr ^O CN e‏ ا at NOx QC o0 -OCN‏ 0۱۱۱۴۵ تر م 
وی ده ته یہ ہ ہے 0 mm) OO OI CIO Ou zt (O00‏ ۳ 00ہ ارہ ن ہہ هو یي ,دص هټېمم سم وہ ۳ ٭ هم 3 — 
0 .9 9 9 9 9 9 9 9 090 9 9 09 99 9 09 9 9 9 0 9 9 099 ق 9 009 9 9 9 9 6 6 9 9 90 9.9 9 9 09690 9 0909 هد ده € 
OR DT CII UD O TINO DOT DIO O CO NI ONO ORIN BR e-01000 0-0 co ON NNS‏ گا ۶ ۷ يہ سم سم سم 

ہی ہے مہو سر Orco‏ 00 0 000 نہ مب NOON) + + + FUN‏ سم سے 

mm) OO r= 


gym 


INDNON DANAN AN DWM AWMS MOVIN ONO NSNDNAN omuwouwaoouyQuYOu OU) Cu Ou OuwuUDOOGOOOO 
€ هد د‎ 6 9 9 00 9 o 6 6 6 6 6 0 6 06 6 6 6 6 6 6 6 6 0 9 o 92 6 9 0 090 6 9 9 9 0 9 9 9 9 9 o یی‎ 6 6 6 6 6 6 oe ه‎ 
NAIM ۴ اټ‎ DO OR OOO ODI خڅ م  لټ )١٤۱ص ۶ - (۳ہ یہ یہ سے سہ رے‎ 00000.6٨ OQ Oe e لیم‎ ۳۲ + +u n O هته مج ن۔‎ OO 

yae pam pae p r— m gm n gam om gm s goma am pem g gamn (NS ENCO CONON CN CNCNCNCN NNN NUN OLN O‏ سم ې 


ÇN]‏ حسم سم 
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THE MARLEX TWO SAMPLE 


APACITIES IN MJ/OK/GM 


~ 
wg 


HEAT 


TARASOV C 


C ACTUAL 


TEMPERATURE 


eee KO TOLO MM OND tm Ore OOM NON TH DO HOM نن‎ — CƏ -O f — CO + —— n CO + C2u n CN -O O — ONO 
e € 9 9 0 می 9 هم‎ 6 560 @  @ o 9 6 60 o 9 6 د‎ o 9 6 6 9 9 0 @? o 0 9 56 هو 6 6 06 9 9 9 9 9 9 ده 6 9 09 یج @ 9 ده ده‎ 
m= ہے ہے ور یر پر‎ | | mM NNO + + ۶ ۱ی +۔‎ ۱۱۱۳(۱۱۸( FT A JF سمسمسریح پت(۶۲۷۱۲۴(۷۷-‎ í | یی ےم‎ QO 
١ 8595/2865 4 9499797 | | | | ٠ ٠ ٠٠ ۱ | | | ١ 48 101 


re NGANG DOO TOO OMA TMM ON SOE ېڼم د سم‎ OME tO EF NNN SOE NOM ORM NWI CO 
ONO O ONO JF ORE OMNI DO -+ ODM NCH MMANDACAYMOM HN TR ATUIOR-B- OMIM AMM r= 
CIN EN CO OUR ېم یم هلص مځ ټم‎ C000 00 O^ e- + OO DM ON ON OM سم یم‎ MDM TNOMOR- OM SMOF-AIOWM 
e e 90 0 099 0 09 90 09 6 6 6 9 9* 0 90 9 ھق‎ 6 6 O 6 6 6 6 6 6 6 6 6 6060606 هه‎ 9 90 6 o 6 6 9 6 6 هه هه‎ ٥ه‎ 9 ٥ هم‎ 
سمسم‎ AIO TID OR DON TUE DN SOA THR 00-900) «9 O roo CQ n0 h- e- zr OO 0) O0 O rA LI OO 0D AICO 

0C 0 OO MIN‏ وه ہد وہ ح ٥مم‏ ہا ضاہ ۱۱۳٣١ ST ۴ ۴ SF‏ ل6 لم لم لم س سم سم سو سم سے 

(NJ O C" —‏ — سے سسس سس مسد 

e 


TF OOM CG Gm (Nr TUN 0 O0 C (C000 zr CJ-O un con P- OF MRR NENAMMOC OG‏ خم تت مده ناواه x X x‏ بې 
AMM DK OOOR-ADNNOAOS HAD ADNOON AAT OFR-DODO OIND—-GOO0DTk-+tOMOo00o0‏ 
هم AMMO R-O OR-OOIN NOOO TS tT ODN ON COOSMONE- DAHDNND HOD OHM OO FOOT‏ 

e © © 0 90 0 9 09 0 0 O O 6 o 9 9 9 09 9 ق می یج 6 ق‎ Oo 9 9 9 9 ده یج 9 09 هه هه‎ 9 6 6 0 6 0 0 9. 09 90 9 @  Ə ə 
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HEAT CAPACITIES IN MJ/OK/GM 
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(AMORPHOUS) SAMPLE 
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